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SECTION  1 
INTRODUCTION 


In  a  continuing  effort  to  provide  support  to  the  Defense  Nuclear 
Agency  (DNA)  infrared  (IR)  program,  particularly  the  field  program,  we  have 
addressed  several  topics,  delineated  in  the  nwork  statement”  to  this  con¬ 
tract,  that  may  be  summarized  as  follows. 

1.  Continue  and  complete  studies  being  done  in  conjunction  with 
AFGL,  Photometries,  and  TIC  to  determine  the  validity  of  our  electron  beam 
deposition  code  by  detailed  comparison  with  EXCEDE: SPECTRAL  optical  data. 

2.  Continue  and  complete  studies  to  determine  the  spatial  dis¬ 
tribution  of  infrared  plasma  radiation  from  electron  beam-irradiated  air, 
and  the  magnitude  of  competing  molecular  band  radiation.  Use  the  results 
to  determine  the  feasibility  of  using  an  EXCEDE-type  experiment  to  check 
the  validity  of  the  theoretical  model  of  plasma  radiation. 

3.  Provide  a  better  estimate  of  the  importance  of  No0  as  an  LWIR 
radiator  in  a  nuclear  environment  by  determining  and  utilizing  better  input 
data  on  its  formation  process  through  the  reaction  N0  (A^E)  +  0^  -*  No0  +  0,  and 
by  attempting  to  determine  the  excitation  rate  of  its  fundamental  bands  by 
sun light  absorpt ion . 

4.  Attempt  to  determine  realistic  rates  for  photon  scatter 
excitation  of  uranium  oxide  and  f'ssion  fragment  oxide  ions  with  the  object 
of  providing  better  estimates  of  the  importance  of  these  species  as  LWIR 
radiators  in  a  nuclear  environment. 
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5.  By  utilizing  prior  MHD  multiburst  calculations,  make  a  pre¬ 
liminary  assessment  of  the  degree  to  which  our  conclusions  regarding  the 
relative  importance  of  LWIR  radiators  in  a  single-burst  environment  are 
modified  in  a  multiburst  environment. 

Item  1,  treated  in  Section  2,  involves  an  attempt  to  define  the 
spatial  dimensions  of  the  anomalously  bright  region  surrounding  the  EXCEDE: 
SPECTRAL  electron  guns,  and  to  determine  if  the  single  particle  deposition 
model  can  possibly  account  for  the  anomaly.  The  usefulness  of  the  EXCEDE 
concept  as  a  diagnostic  tool  for  verifying  our  IR  nuclear-predictive  models 

is  vitally  dependent  on  our  ability  to  understand  the  mechanisms  involved. 

Item  2,  treated  in  Section  3,  involves  a  study,  based  on  ARCTIC 
code  calculations,  of  plasma  radiation  excited  by  a  high-current  beam  of 
3-keV  electrons.  The  main  purpose  of  this  work  is  to  see  if  the  bound- 
bound  component  of  IR  plasma  radiation  in  an  electron  beam  experiment,  such 
as  EXCEDE,  is  likely  to  be  intense  enough,  relative  to  other  competing 

emissions,  to  permit  an  experimental  check  on  the  validity  of  our  theoretical 

models  for  this  radiation. 

Items  3,  4,  and  5,  treated  in  Section  4,  deal  with  questions  that 
arose  in  an  earlier  study  (Reference  1-1)  of  the  relative  importance  of 
various  long  wavelength  infrared  (LWIR)  radiators  in  a  nuclear  environment. 

A  summary  of  our  conclusions  is  presented  in  Section  5. 
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SECTION  2 


ANOMALOUS  BRIGHTENING  IN  THE 
EXCEDE:  SPECTRAL  EXPERIMENT 


PRELIMINARY 

This  section  represents  a  continuation  of  the  work  reported  in 
Section  3  of  Reference  2-1.  It  is  suggested  that  the  reader  become 
acquainted  with  the  details  of  that  section  before  proceeding.  A  brief 
summary  is  as  follows. 

In  order  to  interpret  optical/IR  data  from  electron- impact 

excited  regions  of  the  atmosphere  (e.g.  aurora,  EXCEDE)  in  terms  of  a 

nuclear  environment,  the  tacit  assumption  has  been  that  the  deposition 

mechanism  is  similar  to  that  for  X-ray  photoelectrons,  or  fission  beta 

electrons.  That  is,  the  assumption  has  been  that  the  electrons  deposit 

their  energy  essentially  independently  so  that  collective  effects,  that 

★ 

may  give  rise  to  internal  electric  fields,  are  negligible.  In  such  cases, 
given  the  proper  set  of  cross  sections,  the  electron  energy  spectrum  can  be 
computed,  the  optical/IR  radiation  predicted,  and  our  nuclear-IR  models 
checked  for  accuracy.  If,  on  the  other  hand,  the  electrons  in  a  field 
experiment  should  deposit  their  energy  quite  differently,  because  of  col¬ 
lective  effects,  then  their  energy  spectrum  is  likely  to  differ  from  that 
arising  as  a  result  of  independent  particle  deposition,  as  are  also  the 
relative  efficiencies  for  production  of  various  optical  and  infrared  radi¬ 
ations.  Inferences  drawn  regarding  the  IR  models  used  for  application  to 
the  nuclear  (independent  particle)  case  may  then  be  erroneous. 


*  The  appropriateness  of  this  assumption  for  the  auroral  and  nuclear 
environments  is  treated  in  a  later  subsection. 
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In  order  to  check  the  applicability  of  the  independent  particle 
model  to  EXCEDE  experiments,  we  used  a  computer  code,  that  incorporates 
such  a  model,  to  calculate  the  energy  deposition  rate  surrounding  electron 
beams  in  the  atmosphere.  Applying  a  fluorescence  efficiency  factor  of 

-  3  +  o  * 

4.5 x  10  for  N0  1  Neg.  (3914  A)  emission,  we  then  compared  the  results 
with  photographic  data  from  the  PRECEDE  II  and  EXCEDE : SPECTRAL  events. 

The  results  of  those  comparisons  were  reported  in  Reference  2-1. 

Briefly,  the  PRECEDE  II  comparisons  were  somewhat  inconclusive 
because  absolute  intensities  are  not  available  from  the  field  data.  However, 
the  EXCEDE : SPECTRAL  data  are  greatly  at  odds  with  the  code  calculations  and 
seem  to  imply  that  the  independent  particle  model  is  woefully  inadequate  to 
account  for  the  observations,  especially  in  a  region  within  a  few  tens  of 

o 

meters  of  the  electron  accelerator.  In  this  region,  the  measured  5914  A 
brightness,  which  is  proportional  to  the  energy  deposition  rate,  is  up  to 
two  orders -of -magnitude  larger  than  the  calculated  values. 

Several  questions  immediately  arise.  For  example:  What  is  the 
mechanism  that  produces  the  anomalous  brightening?  Does  it  involve  collec¬ 
tive  effects  that  may  render  interpretation  of  optical/IR  data  questionable? 
Over  what  region  does  the  anomaly  pertain?  Is  there  a  region  sufficiently 

far  from  the  accelerator  where  the  deposition  behaves  ’’normally?"  Is  the 
+  ° 

First  Negative  (3914  A)  fluorescence  efficiency  inferred  from  auroral 
data,  and  measured  under  low-current  conditions  in  the  laboratory,  inappli¬ 
cable  under  the  high-current  conditions  of  EXCEDE?  Answers  to  most  of  these 
questions  are  not  known  for  certain  at  the  present  time,  although  we  will 
attempt  to  address  some  of  them  in  this  section. 


*  This  factor  is  consistent  with  auroral  data  and  with  many  laboratory 
experiments  involving  low-current  (usually  in  the  tens  of  ya  range) 
electron  beams  (Reference  2-2). 
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Possible  mechanisms  that  have  been  suggested  include  various  dis¬ 
charge  phenomena,  particularly  the  so-called  ’’beam  plasma  discharge"  (BPD) 
observed  in  laboratory  experiments  (References  2-3  to  2-6),  return  currents 
induced  by  vehicle  charging,  and  payload  outgassing.  The  latter  mechanism, 
however,  probably  plays  a  relatively  minor  role  since  the  anomaly  is  present 
even  at  altitudes  as  low  as  80  km  where  the  ambient  density  is  believed  to 
be  considerably  greater  than  that  from  outgassing.  At  present,  some  dis¬ 
charge  phenomenon  appears  to  be  favored. 

If,  in  fact,  a  BPD  is  actually  responsible,  it  is  noteworthy  that 
vacuum  tank  experiments  (Reference  2-7)  have  shown  that  the  electron  energy 
spectrum  is  then  dramatically  changed,  with  electron  energies  greater  than 
the  incident  beam  energy  present  along  with  a  greatly  enhanced  low  energy 
tail.  As  mentioned  above,  one  would  then  have  to  view  any  extrapolation  of 
optical/IR  data  therefrom  to  the  nuclear  case  with  considerable  suspicion. 
Alternatively,  validation  of  our  optical/IR  models  could  still  be  accom¬ 
plished  provided  we  knew  the  electron  energy  spectrum  throughout  the  dosed 
region.  However,  this  would  require  detailed,  and  probably  difficult, 
measurements  which  should  be  avoided  if  at  all  possible. 

For  interpretation  of  the  EXCEDE: SPECTRAL  data,  and  particularly 
for  use  in  planning  any  future  EXCEDE  experiment,  we  would  especially  like 
to  know:  (1)  the  extent  of  the  anomalous  region,  and  (2)  whether  or  not 
the  electron  energy  spectrum  in  this  region  differs  significantly  from  that 
predicted  on  the  basis  of  independent  particle  deposition.  Closely  related 
to  Item  2  is  the  question  whether,  under  high-current  conditions,  the 

o 

effective  fluorescence  efficiency  for  3914  A  emission  can  be  raised  to  the 
point  where  the  observations  can  still  be  accounted  for  on  the  basis  of 
independent  particle  theory.  This  section  addresses  these  two  topics. 
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FURTHER  COMPARISONS  WITH  FIELD  DATA 


As  in  Reference  2-1,  Section  3,  the  data  and  calculations  dis¬ 
cussed  here  all  pertain  to  an  electron  pulse  initiated  at  an  altitude  of 
about  123  km  on  rocket  ascent  during  the  EXCEDE: SPECTRAL  event. 

On  Size  of  Anomalous  Region 

The  best  hope  of  estimating  the  size  of  the  enhanced  region  comes 
from  a  study  of  the  ground-based  photographic  data  taken  by  Technology 
International  Corporation  (TIC).  Figure  3-9  of  Reference  2-1  shows  a  com¬ 
parison  between  the  relative  radiance  values  measured  by  TIC,  as  a  function 

o 

of  distance  along  the  beam,  and  our  code  calculations  for  3914  A.  In  that 
figure  we  had  arbitrarily  normalized  the  TIC  data  to  correspond  to  the 
calculated  curve  at  a  point  3  km  along  the  beam  from  the  accelerator. 
Actually,  however,  the  figure  is  not  particularly  meaningful  for  two 
reasons.  First,  because  an  absolute  intensity  calibration  for  the  data  was 
not  available  at  that  time  and,  second,  because  the  spatial  resolution  of 
the  data  is  estimated  to  be  in  the  range  from  about  20  to  50  meters  (Refer¬ 
ence  2-8),  whereas  the  theoretical  curve  corresponds  to  infinitely  sharp 
(0  meter)  spatial  resolution. 

A  somewhat  more  meaningful  comparison  is  shown  below  in  Figure  2-1. 
Here,  the  preliminary  values  of  absolute  intensity  measured  by  TIC  (Refer¬ 
ence  2-8  are  compared  with  four  different  theoretical  curves  corresponding 
to  spatial  resolutions  of  0,  10,  20,  and  30  meters.  Since  the  resolution 
of  the  camera  data  is  believed  to  exceed  20  meters,  the  lowest  of  the  (solid) 
theoretical  curves  is  probably  the  most  appropriate  for  comparison  purposes. 
Nevertheless,  there  is  still  a  certain  degree  of  uncertainty  in  making  this 

comparison,  partly  because  of  the  resolution  question,  and  partly  because 

0  + 

the  theoretical  curves  refer  to  radiance  in  just  the  3914  A  (N^  1  Neg.) 

band  whereas  the  photographic  data  include  other  wavelengths  as  well. 
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Figure  2-1.  EXCEDE: SPECTRAL  radiance  (3914A)  normal  to  beam  axis  at 
123  km;  calculations  and  photographic  data  compared. 


Thus,  to  compare  the  data  and  calculations  on  an  equal  basis,  the  data 
curve  (dashed)  should  be  lowered  by  some  factor,  estimated  to  be  about  1.5. 

A  still  more  meaningful  comparison  can  be  made  on  the  basis  of 
the  crosswise  integrated  intensity  as  a  function  of  distance  along  the 
beam.  This  has  the  advantage  of  removing  the  dependence  of  the  results  on 
the  spatial  resolution  of  the  observing  equipment.  The  crosswise  integrated 
intensity,  T(Z),  can  be  defined  as  the  integrated  intensity  along  a  line 
normal  to  the  beam  at  a  distance  Z  from  the  accelerator.  That  is,  if 
I(Z,R)  is  the  intensity  normal  to  the  beam  at  a  distance  Z  along  the  axis 
from  the  accelerator,  and  at  a  radial  distance  R  from  the  axis,  then 

00 

T(Z)  = J  I(Z,R)dR  .  (2-1) 

-00 

Thus,  given  a  set  of  radial  intensity  distributions,  I(Z,R),  at  various 
distances  along  the  beam  from  the  accelerator,  the  area  under  each  gives 
the  crosswise  integrated  intensity. 

Figure  2-2  shows  the  calculated  (lower  solid  curve)  and  observed 
(dashed  curve)  radial  intensity  distributions  at  a  distance  of  4  meters 
along  the  beam  from  the  gun.  The  data  were  taken  by  Photometries  (Refer- 
ence  2-8)  using  a  rocket-borne  black  and  white  camera  whose  filter  band 
extends  from  about  3600  X  to  4700  A.  However,  since  the  calculated  curve 

o 

is  for  the  3914  A  band  only,  the  dashed  curve  should  be  lowered  by  some 
factor,  estimated  to  be  about  1.5  (Reference  (2-10).  Even  so,  it  is  clear 
that  the  area  under  the  dashed  curve  (crosswise  integrated  intensity)  is 
considerably  greater  than  that  under  the  theoretical  curve.  The  upper 
solid  curve  in  Figure  2-2  is  the  calculated  intensity  distribution 
increased  by  a  factor-  of  10  and  shifted  to  the  left  to  bring  its  peak  into 
line  with  the  camera  data.  More  will  be  said  about  that  in  a  later  sub¬ 
section  . 
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Figure  2-2.  Radial  intensity  distribution  at  +  4.0  meters  from  EXCEDE: 
SPECTRAL  gun  (123-km  altitude). 


Figures  2-3  to  2-5  are  examples  of  similar  comparisons  between 
the  calculated  and  observed  radial  intensity  distributions  at  distances 
along  the  beam  from  the  gun  of  44,  324,  and  1450  meters,  respectively. 

In  these  figures  the  data  are  from  the  ground-based  photographic  image 
recorded  by  TIC  (Reference  2-11).  Again,  as  mentioned  above,  since  the 
data  includes  all  wavelengths  for  which  the  camera  film  was  sensitive,  the 
dashed  curves  should  be  lowered  by  some  factor  (*^1.5).  Nevertheless,  a 
visual  inspection  of  Figures  2-3  to  2-5  shows  that  the  area  under  the 
observed  intensity  distribution  (crosswise  integrated  intensity)  is  greater 
than  than  under  the  theoretical  distribution  at  44  meters  from  the  gun 
(Figure  2-3),  that  the  reverse  is  true  at  1450  meters  (Figure  2-5),  and 
that  the  two  areas  are  equal  at  some  intermediate  distance. 

The  crosswise  integrated  intensity,  defined  by  Equation  2-1,  can 
be  defined,  alternatively,  (for  an  optically  thin  medium)  in  terms  of  the 

O  • 

volume  emission  rate  at  3914  A,  by  the  expression 

00 

T3914(Z)  =  j27rRe3914(Z,R)dR  (2-2) 

o 

or,  in  terms  of  the  volume  deposition  rate,  e(Z,R),  in  normalized  units 
of  eV  gm  1  (See  Reference  2-1),  by  the  expression 

oo 

^3914^  =  J* Re(Z,R)dR  (ergs  cm  *  sec  *  ster  ^)  .  (2-3) 

o 

In  Equation  2-3,  I  is  the  beam  current  in  amperes,  p  is  the  local  atmo- 

-3 

spheric  density  (gm  cm  )  and  the  radial  distance,  R,  is  measured  in  meters. 

Equations  2-2  and  2-3  were  used,  together  with  values  for  e(Z,R) 
computed  by  our  deposition  code  (see  Reference  2-1,  Figure  3-3),  to  calcu¬ 
late  the  crosswise  integrated  intensity  as  a  function  of  distance  along  the 
beam.  The  results  are  shown  in  Figures  2-6  and  2-7.  Figure  2-6  gives  the 
results  in  normalized  units  for  a  single  3-keV  electron  whereas  Figure  2-7 
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Figure  2-3.  Radial  intensity  distribution  at  +  44  meters  from  EXCEDE: 
SPECTRAL  gun  (123-km  altitude). 
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Figure  2-4.  Radial  intensity  distribution  at  324  meters  from  EXCEDE: 
SPECTRAL  gun  (123-km  altitude). 
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Figure  2-5.  Radial  intensity  distribution  at  +  1450  meters  from  EXCEDE: 
SPECTRAL  gun  (123-km  altitude). 
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Figure  2-6.  Calculated  crosswise  integrated  intensity,  T(Z),  in  magnetic  field 
of  0.5  gauss  in  normalized  units  for  unit  number  current  exiting  gun 
in  30°  cone  with  axis  parallel  to  §  (E  =3  keV). 


+  CQ 
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Figure  2-7.  Calculated  crosswise  integrated  intensity  (3914X)  at  123  km  for 
7  amp,  3  keV  beam  exiting  gun  in  30°  cone  with  axis  parallel  to 


shows  the  results  for  the  particular  case  of  interest  here,  namely  a  7  amp 
beam  at  123-km  altitude.  In  both  figures  the  electron  gun  is  located  at 
the  origin,  and  emission  for  negative  values  of  the  distance  arises  from 
backscattered  electrons. 

Figure  2-8  compares  the  calculated  crosswise  integrated  intensity 
(solid  curve)  with  that  observed  (dashed  curve).  The  observed  values  were 
obtained  by  finding  the  areas  under  the  radial  intensity  distributions  pro¬ 
vided  by  TIC.  For  this  purpose,  we  arbitrarily  extrapolated  the  wings  of 

* 

the  distributions  to  zero  intensity.  A  value  for  T(Z),  inferred  from  the 
Photometries  data  (Figure  2-2)  for  a  distance  of  4  meters  forward  of  the 
gun,  is  also  shown  in  Figure  2-8. 

A  comparison  of  the  twro  curves  shown  in  Figure  2-8  provides  some 
measure  of  the  extent  of  the  enhanced  region.  From  the  points  of  inter¬ 
section  of  the  dashed  and  solid  curves,  we  infer  that  the  enhancement 
extends  from  about  280  m  behind  the  gun  to  about  800  m  in  front  of  it. 

These  values,  however,  must  be  changed  somewhat  if  we  allow  for  the  fact 
that  the  dashed  curve  should  be  lowered  by  some  factor,  as  mentioned  above, 

o 

because  it  represents  emission  over  a  broader  band  than  just  the  3914  A 
band  represented  by  the  theoretical  curve.  For  example,  if  we  assume  the 
factor  is  1.5,  the  enhancement  then  extends  from  about  200  m  behind  to 
about  550  m  in  front  of  the  gun.  If  the  factor  is  2.0,  it  extends  from 
about  130  m  behind  to  about  350  m  in  front. 

Further  considerations,  introduced  in  the  following  subsection, 
lead  to  the  conclusion  that  the  extent  of  the  enhanced  region,  as  inferred 
from  the  TIC  data,  is  even  less  precisely  defined  than  the  foregoing  numbers 
suggest . 


*  The  implication  of  these  extrapolations  to  the  accuracy  of  the  results 
is  discussed  in  the  next  subsection. 
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Tigure  2-8.  Crosswise  integrated  intensity  at  123-km  altitude 
(EXCEDE:SPECTRAL) . 
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On  Energy  Conservation  and  Accuracy  of  Above  Results 


In  Figure  2-6  we  showed  the  calculated  crosswise  integrated  in¬ 
tensity,  in  normalized  units,  for  a  current  of  1  (3  keV)  electron  sec  ^ 
exiting  the  gun  in  a  30°  cone  with  axis  parallel  to  the  earth1 s  magnetic 
field.  Hie  area  under  this  curve  (extended  from  -0oto+ool  equals  3  keV. 
However,  for  a  magnetic  dip  angle  of  77°  (Poker  Flat,  Alaska),  the  column 

density  of  air  along  the  earth !s  field  above  an  altitude  of  123  km  is  about 
-6  -  2 

23  x  10  gm  cm  .  Thus,  in  Figure  2-6,  the  energy  represented  by  that 

_  6  _  o 

portion  of  the  curve  for  which  Z>23  xio  gm  cm  would  escape  from  the 
atmosphere.  From  the  respective  areas  involved,  we  find,  on  the  basis  of 
the  independent  particle  model,  that  about  8  percent  of  the  energy  should 
escape  the  atmosphere. 

Turning  now  to  the  actual  EXCEDF: SPECTRAL  observat ions ,  it  is 

of  interest  to  see  how  the  power,  inferred  from  the  area  under  the  dashed 

curve  in  Figure  2-8,  compares  with  the  total  power  in  the  EXCEDE  beam. 

For  this  purpose,  we  have  arbitrarily  extrapolated  the  dashed  curve  to  zero 

in  both  wings,  measured  its  area,  and  multiplied  by  4tt.  The  result  is 
8  - 1 

1.43  x  10  ergs  sec  in  the  photographic  band.  If  we  assume  that  the  photo¬ 
graphic  film  was  sensitive  to  50  percent  more  energy  than  represented  by 

o  o 

just  the  3914  A  band,  and  that  the  fluorescence  efficiency  of  the  3914  A 
-3 

band  is  4.5  x  10  ,  then  the  total  power  represented  by  the  TIC  data,  is 

O  *7  If)  _  1 

PTIC  =  1.43  x  10/(1 . 5  x  4.5  x  10~j  =  2.1X101  (ergs  sec  *(2-4) 

But  the  power  in  the  electron  beam  is 

^beam  =  (electrons/sec“aiT1P)  x  (amps)  x  (eV/electron)  x  (ergs/eV) 

=  (6. 25  x  1018)  x  (7)  x  (3  x  103)  x  (1.6  x  10"12) 

=  2.  lx  10^  (ergs  sec  L  .  (2-5) 
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Therefore,  the  TIC  data  seem  to  account  for  only  10  percent  of  the  beam 
energy.  Where  has  the  rest  of  it  gone? 

For  one  thing  part  of  it,  but  probably  no  more  than  about  8  percent 

★ 

as  indicated  above,  escaped  the  atmosphere.  However,  most  of  the  missing 

energy,  we  believe,  is  lost  in  the  film  fog.  That  is,  it  was  deposited  in 
large  volumes  of  air  outside  the  periphery  of  the  visible  image  such  that 
the  irradiancy  at  the  camera,  from  those  regions,  was  below  the  film 
threshold. 


This  raises  a  question  concerning  the  accuracy  of  our  determination 
of  the  crosswise  integrated  intensity  from  the  measured  radial  intensity 
distributions.  This  is  because  the  wings  of  the  measured  curves  (see 
dashed  curves  in  Figures  2-3  to  2-5)  are  below  film  threshold  and  we  can 
only  guess  how  they  fall  off  at  larger  radial  distances.  In  obtaining  the 
dashed  curve  in  Figure  2-8,  we  used  a  linear  extrapolation  of  the  wings  to 
zero  intensity  for  each  of  the  measured  radial  distributions.  More  likely 
the  intensity  in  the  wings  falls  off  less  rapidly,  as  it  does  in  our  indepen¬ 
dent  particle  model,  with  the  result  that  each  point  in  the  dashed  curve  of 
Figure  2-8  should  be  raised  by  some  unknown  factor.  Similarly,  for  the 
purpose  of  measuring  the  area  under  the  dashed  curve  in  Figure  2-8,  we 
linearly  extrapolated  the  two  ends  to  zero.  They,  too,  probably  fall  off 
less  rapidly  with  distance. 

The  conclusion  then  is,  that  in  order  to  conserve  energy,  the 
dashed  curve  in  Figure  2-8  needs  to  be  raised  by  some  unknown  amount,  and 
tails  added  that  extend  to  great  distances  in  both  forward  and  backward 
directions.  It  is  not  possible,  however,  from  the  photographic  data,  to 


*  One  would  expect  that  since  an  abnormally  large  fraction  of  the  energy 
is  deposited  near  the  electron  gun,  a  smaller  than  normal  fraction 
would  be  available  to  escape. 
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quantify  these  changes.  All  we  can  say  regarding  the  dashed  curve  in 
Figure  2-8  is  that  it  should  probably  be  raised.  With  reference  to  the 
preceding  subsection,  this  means  that  we  cannot  be  at  all  precise  about 
the  spatial  dimensions  of  the  anomalous  region  surrounding  the  accelerator, 
other  than  to  say  that  it  appears  to  extend  for  several  hundred  meters  along 
the  magnetic  field. 

HOT  PLASMA  EXCITATION  OF  3914  A 

In  the  preliminary  discussion  of  this  section,  we  mentioned 
several  mechanisms  that  have  been  suggested  as  possible  candidates  to 
explain  the  anomalous  enhancement  of  emission  observed  in  the  EXCEDE : SPECTRAL 
experiment.  Favored  among  these  mechanisms  is  some  discharge  phenomenon, 
possibly  the  ’’beam  plasma  discharge"  (BPD)  that  has  been  studied  rather 
extensively  in  the  laboratory.  It  was  also  mentioned  that  the  latter 
mechanism  results  in  an  electron  energy  spectrum  that  deviates  markedly 
from  what  one  expects  from  electron  energy  deposition  in  the  aurora,  or 
following  a  nuclear  detonation,  for  which  the  independent  particle  model 
should  be  a  very  good  approximation.  Since  the  efficiency  for  production 
of  certain  optical  and  IR  bands  may  be  sensitive  to  the  electron  energy 
spectrum,  use  of  optical/IR  data  from  discharge-excited  air  to  validate 
our  nuclear-IR  models  may  be  hazardous.  Of  course,  given  the  electron 
energy  spectrum  we  can  use  codes,  such  as  ARCTIC,  to  predict  the  optical/IR 
radiation  for  comparison  with  the  data,  and  thereby  check  the  models.  But, 
"given  the  electron  energy  spectrum"  means  measuring  it  throughout  the 
dosed  region,  which  is  a  difficult  and  probably  expensive  thing  to  do. 
Therefore,  insofar  as  EXCEDE  is  to  be  considered  a  viable  tool  for  helping 
to  validate  our  nuclear-IR  codes,  it  would  be  of  great  benefit  if  it  could 
be  shown  that  the  observed  enhancement  can  be  explained  on  the  basis  of 
independent  particle  deposition.  Partly  toward  that  end,  the  work  in  this 
subsection  was  undertaken. 
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The  only  way  we  can  think  of  whereby  the  independent  particle 
model  might  be  able  to  account  for  the  anomalous  brightening  is  if  the 

o 

effective  fluorescence  efficiency  for  3914  A  radiation  is  greatly  increased 
by  the  presence  of  plasma  electrons  that  have  been  heated  by  collisions  with 
beam  electrons  to  the  point  where  they  can  produce  further  ionization  and 

o 

electronic  excitation.  The  fluorescence  efficiency  for  3914  A  that  we  have 
used  in  our  calculations,  reported  in  Reference  2-1  and  earlier  in  this 
section,  is  based  on  the  results  of  numerous  laboratory  measurements  (see, 
for  example,  the  references  given  in  Reference  2-12)  as  well  as  inferences 
from  auroral  data.  However,  as  far  as  we  can  determine,  in  all  of  those 
experiments  the  current  of  electrons  involved  was  too  small  to  create  a 
plasma  hot  enough  to  produce  significantly  more  ionization.  Therefore,  to 
our  knowledge,  no  experiment  or  calculation  has  previously  been  done  to 
determine  how  plasma  heating  by  a  high-current  (amperes)  electron  beam  in 
air  might  act  to  modify  the  effective  fluorescence  efficiency. 

To  investigate  this,  we  have  made  use  of  the  ARCTIC  code.  This 
code  was  originally  constructed  by  MRC,  partly  under  DNA  auspices,  for  use 
in  analyzing  and  interpreting  optical/IR  data  from  the  auroral  ICECAP 
experiments.  Details  of  the  code,  for  auroral  use,  can  be  found  in  Refer¬ 
ences  2-13  through  2-16,  and  a  very  brief  summary  is  given  in  the  following 
subsection.  For  use  under  high-current  conditions,  however,  we  found  it 
necessary  to  make  extensive  changes  in  the  code.  A  detailed  description 
of  these  changes  is  presented  in  Appendix  A,  and  a  brief  summary  is  given 
immediately  following  the  next  subsection.  The  changes,  that  include  the 
introduction  of  additional  cooling  mechanisms,  were  necessary  not  only  for 
the  investigations  of  this  section  but  also  for  the  study  of  plasma  radia¬ 
tion  in  EXCEDE  experiments,  described  in  Section  3. 
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Summary  of  Auroral  Code  "ARCTIC 


ARCTIC  is  a  general  purpose  atmospheric  electron  deposition  and 
chemistry  code  originally  used  to  study  the  effects  of  a  flux  of  auroral 
electrons  upon  the  atmosphere.  It  calculates  the  electron  energy  deposition 
and  partition,  the  subsequent  atmospheric  chemistry,  optical  and  IR  emissions, 
and  electron  heating  effects. 

Specifically,  we  can  input  a  spectral  flux  of  primary  electrons 
with  a  given  pitch-angle  distribution  and  a  specified  magnetic  dip  angle. 

The  electrons  then  interact  with  the  atmosphere,  consisting  of  N0,  0o,  0, 

CO^,  N,  and  an  ambient  electron  gas.  The  code  computes  a  number  of  items, 
as  a  function  of  altitude,  including  the  secondary  electron  production  rate 
spectrum,  secondary  flux  spectrum,  and  the  volume  production  rates  of 
excited  states  by  both  the  primaries  and  by  all  generations  of  secondary 
electrons . 


The  production  rates  of  allowed  states,  i.e.,  the  non-metastables, 
are  then  used,  with  allowance  for  branching  and  cascade,  to  compute  directly 
the  volume  emission  rates  in  various  allowed  lines  and  bands.  The  production 
rates  of  metastable  states  and  the  electron  gas  heating  rate,  computed  in  the 
deposition  routine,  are  then  fed  into  the  chemistry  subroutine  where  the 
species  concentrations  and  electron  and  vibrational  temperatures  are 
calculated.  These  quantities  are  then  used  to  compute  the  volume  emission 
rates  of  radiation — both  optical  and  IR  —  resulting  from  forbidden  transi¬ 
tions  and  also  from  vibraluminescent  and  chemiluminescent  reactions.  These 
emission  rates  are  output  as  functions  of  altitude  and  time. 

In  particular,  the  code  provides  for  the  excitation  by  the  bom¬ 
barding  flux  of  210  states  (including  No(v=0->7))  among  the  species  N0,  N2+ , 

N,  N’+ ,  02>  0„+,  0,  0+,  CO,, ,  and  C0o+,  while  the  chemistry  subroutine 
treats  over  200  reactions  among  31  chemical  species  whose  concentrations 
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are  determined  by  integration  of  the  rate  equations  using  a  combination 
Runge-Kutta  and  quasi-analyt ic  approach.  For  determination  of  the  electron 
temperature,  the  code  includes  direct  heating  by  the  bombarding  electron 
flux  and  cooling  by  excitation  of  the  metastable  states  of  0,  0+,  and  N, 

rotational  and  vibrational  excitation  of  N  and  0  ,  fine  structure  excita- 

3  l  l 

tion  of  0(  P) ,  and  elastic  collisions  with  0^,  0,  and  positive  ions. 

Volume  emission  rates  output  from  the  code  include  the  following 
bands.  From  direct  excitation  of  allowed  states:  the  N0  +  First  Negative 
and  Meinel  bands,  the  Second  Positive,  0?  First  Negative,  and  certain 
IR  bands  of  CO^.  From  the  species  carried  in  the  chemistry  code:  NO  and 
N0+  chemiluminescence  in  the  fundamental  and  first  overtone  bands,  C0o 
vibraluminescence  in  various  IR  bands,  forbidden  line  radiation  from  0,  0+, 
N,  and  N+,  and  bands  of  the  0o  atmospheric  and  infrared  atmospheric  systems. 
Full  provision  for  the  effects  of  quenching  are  included. 

Summary  of  ARCTIC  Code  Changes  Required 
Under  High-Current  Conditions 

Under  high-current  conditions,  we  expect  the  electron  temperature 
to  attain  values  much  higher  than  those  in  the  aurora.  It  then  becomes 
necessary  to  allow  for  cooling  of  the  plasma  electrons  by  other  processes 
including  excitation  of  some  of  the  higher-lying  electronic  states,  and 
even  ionization  of  the  major  species.  Furthermore,  while  attempting  to  run 
the  code  under  EXCEDE-bombardment  conditions,  we  found  that  certain  cooling 
rate  expressions,  valid  under  auroral  conditions,  were  invalid  at  higher 
electron  temperatures.  This  necessitated,  in  those  cases,  drastic  revision 
of  the  rates.  The  details  of  these  and  other  changes  that  have  been  made 
are  described  in  Appendix  A.  A  brief  summary  is  as  follows. 
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Additions  to  ARCTIC  Code 


•  Electron  cooling  by  electronic  excitation  and  ionization  of 
24  additional  states  of  ,  0o,  0. 

o 

•  Emission  of  3914  A  photons  through  the  processes: 

•  ep  +  N2  CX1^)  -  2  ep  +  N'2+(R2I) 

+  2  0 
N0  (X  E)  +  hv (3914  A) 

•  ep  +  N2+(X2I)  -  ep  +  N2+(R2X) 

->  N2+(X2E])  +  hv (3914  A) 

where  e^  represents  a  plasma  electron. 

•  Ionization  of  (v>0) 

•  Ionization  of  O^fa^A)  and  0n(b^Z) 

•  Radiative -col lisional  recombination 

•  Changes  in  ARCTIC  Code  Models 

•  Electron  cooling  by  elastic  scattering  from  N  ,  0^,  0 

•  Electron  cooling  by  vibrational  excitation  of  N0 ,  0^. 

Code  Calculations  and  Results 

After  making  the  required  changes  and  additions  to  the  ARCTIC 

code,  we  exercised  it  under  input  conditions  approximately  simulating  those 

in  the  EXCEDE: SPECTRAL  experiment  at  an  altitude  of  123  km.  The  volume 

12  -  3  - 1 

deposition  rate  was  taken  to  be  1.48x10  eV  cm  '  sec  ,  which  equals  the 
maximum  value  calculated  (using  the  independent  particle  model)  for  a  7  amp, 
3  kV  beam  at  123-km  altitude  (see  Reference  2-1,  Figure  3-3).  The  code  was 
run,  under  full  bombardment  conditions,  for  times  to  10  sec. 
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The  code  output  is  voluminous.  It  includes  production  rates  for 
over  200  states  as  well  as  the  time  dependence  of  more  than  30  species  con¬ 
centrations,  the  electron  temperature,  and  numerous  optical/IR  volume 
emission  rates.  Figure  2-9  shows  only  that  part  of  the  output  that  is 
especially  relevant  to  this  discussion.  Included  is  the  electron  temperature 

o 

and  four  different  contributions  to  the  volume  emission  rate  of  3914  A  for 
-4 

times  from  10  sec  to  10  sec.  These  contributions,  shown  by  curves  2,  3, 

4,  and  5  in  Figure  2-9,  arise  from,  respectively,  (a)  direct  bombardment  of 

NL  by  beam  electrons,  (b)  direct  bombardment  of  N  +  by  beam  electrons,  (c) 

^  ^  + 
excitation  of  N9  by  hot  plasma  electrons,  and  (d)  excitation  of  N~  by  hot 

l  g  1 

plasma  electrons.  The  dashed  line  1,  corresponding  to  2.07*  10  photons 

-3-1  .  0 

cm  sec  ,  is  based  on  a  fluorescence  efficiency  for  the  3914  A  band  of 

-  s 

4.5x  10  ,  which  is  the  value  used  in  our  comparisons  with  the  EXCEDE : SPECTRAL 

photographic  data. 

°  + 

Process  (b) ,  excitation  of  3914  A  by  bombardment  of  Nn  by  beam 

electrons,  is  not  actually  in  the  ARCTIC  code.  We  have  estimated  its  effect 

(curve  3  in  Figure  2-9)  by  noting  that  the  cross  section  for  the  process  is 

4  2 

approximately  8.9  times  that  for  the  process  +  S)-*e^^+N(  P)  > 

which  is  in  the  code.  Thus,  the  volume  production  rate,  P(N^+(R),  of  the 
parent  state  of  the  First  Negative  system,  by  electron-beam  bombard¬ 

ment,  can  be  written  as 


P  (N  (B)  ),  M  + 

K  2  ^beam  on 


^  8.9  P(N(2P))  [N2+] 


(2-6) 


2  .24 

where  p(N(  P))  is  the  production  rate  of  N(  P)  atoms  per  N(  S)  atom.  But, 

from  Reference  2-12,  the  volume  emission  rate,  ^ 3 9 24 »  the  (0>  0)  band 

(3914  A)  of  the  Nn+  First  Negative  system  is  given  by 


^3914  “  0,65  P(N2+(B)) 


(2-7) 
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(Photons  cni  sec 


10 


10  c  10' 
Bombardment  Time  (sec) 


Figure  2-9.  Calculated  3914^  volume  emission  rate  at  123-km 

altitude  for  an  energy  deposition  rate  of  1.48x10 u 
eV  cm-3  sec-1  by  a  3  keV  electron  beam. 


32 


Electron  Temperature  (°K) 


Therefore, 


0^3914*^  ^  1.6  [N?  +  ]  (photons  cm  3  sec  *)  (2-8) 

beam  on 

2 

where,  for  p(N(  P)),  we  have  used  the  ARCTIC  code  value  of  0.276. 

A  study  of  the  curves  in  Figure  2-9  reveals  the  following.  Until 

a  time  of  0.6  sec,  the  dominant  emission  mechanism  is  the  direct  ionization 

and  excitation  of  N0  by  the  beam  electrons.  For  times  between  0.6  sec  and 

about  3.5  sec,  however,  the  electron  temperature,  the  electron  density,  and 

the  N?  concentration  are  sufficiently  large  that  plasma-electron  excitation 

of  N0+  is  dominant.  For  times  to  about  0.3  sec,  the  second  most  important 

mechanism  (based  on  our  approximate  treatment  of  it  as  outlined  above)  is 

excitation  of  N0+  ions  by  the  beam  electrons.  At  no  time  are  the  electron 

temperature  and  electron  density  simultaneously  high  enough  for  plasma- 

electron  ionization  and  excitation  of  N?  to  be  important.  The  total  emission 

rate  (curve  6)  exceeds  the  nominal  value  (curve  1)  by  a  factor  of  about  2.5 

for  bombardment  times  between  about  1  and  2  sec,  but  after  about  3.6  sec 

it  is  less  than  the  nominal  value  because,  by  this  time,  the  N0  concentra- 
*  ^ 
tion  has  decreased  by  dissociation. 

o 

Although,  as  mentioned  above,  the  maximum  increase  in  the  3914  A 
emission  rate  beyond  the  nominal  value  is  a  factor  of  about  2.5,  the  in¬ 
crease  during  the  actual  bombardment  time  of  an  air  parcel  in  the  EXCEDE: 
SPECTRAL  experiment  (which  is  certainly  less  than  about  0.1  sec  because  of 
rocket  motion)  is  seen  from  Figure  2-9  to  be  no  more  than  about  50  percent. 
This  is  a  far  cry  from  the  2  orders-of-magnitude  or  more  increase  required 
to  account  for  the  anomalous  enhancement. 


*  See  Figure  3-1  of  Section  3. 
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The  conclusion  then  is  that,  even  for  high-current  beams,  plasma 
heating  by  electrons  that  behave  like  independent  particles  is  insufficient 

o 

to  raise  the  effective  fluorescence  efficiency  at  3914  A  to  the  point 
where  it  can  account  for  the  observed  EXCEDE : SPECTRAL  data. 

SOME  BEAM  PLASMA  DISCHARGE  CONSIDERATIONS 

As  mentioned  earlier  in  this  section,  a  likely  candidate  to 
explain  the  observed  enhancement  in  the  EXCEDE : SPECTRAL  experiment  is  the 
beam  plasma  discharge  (BPD)  phenomenon  that  has  been  studied  in  the  labora¬ 
tory  under  both  pulse  (References  2-3,  2-4)  and  DC  (References  2-5  and  2-6) 
conditions,  and  may  have  been  observed  in  space  during  the  Norwegian  Polar  5 
experiment  (Reference  2-17).  In  this  subsection  we  briefly  summarize  some 
of  the  criteria  that  appear  necessary  for  onset  of  a  BPD,  together  with  the 
observed  optical  consequences,  and  then  attempt  to  determine  the  applica¬ 
bility  of  the  phenomenon  to  nuclear  burst,  auroral,  and  EXCEDE  environments. 

Laboratory  Observations  and  Inferences 

The  large  vacuum  tank  experiments  of  Bernstein  et  al.  can  be 

interpreted  (Reference  2-6)  as  defining  a  critical  current,  I  ,  for  onset 

-5  C 

of  a  BPD,  for  pressures  P<  2x  10  torr,  by  the  expression 

-10  v3/2 

I  ~  10  ^ —  Camps)  .  (2-9) 

B  PL 

Here,  V  is  the  beam  energy  (volts),  B  is  the  magnetic  field  (gauss),  and 

L  is  the  system  length  (meters)  which,  in  the  Bernstein  experiments  was  20  m. 

Typically,  for  a  2  keV  beam,  I  was  found  to  have  a  value  of  about  40  ma 

which  implies,  for  an  effective  beam  diameter  of  1  meter,  a  critical  current 

-6  -2  * 

density  of  about  5*  10  amps  cm 

*  This  value  for  the  critical  current  density  is  also  consistent  with  the 
critical  beam  electron  density  shown  in  Equation  2-14. 
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Alternatively,  the  Bernstein  experiments  (References  2-5,  2-6), 
run  under  DC-current  conditions,  and  also  the  Getty  and  Smullin  experiments 
(Reference  2-3),  run  under  pul sed-current  conditions,  indicate  that  a  neces¬ 
sary  condition  for  BPD  ignition  occurs  when  the  beam  current  produces  a 
critical  plasma  density,  N  ,  such  that  the  plasma  frequency,  u)^,  and  the 

More  recent  experiments 


cyclotron  frequency,  go  ,  are  approximately  equal 

(Reference  2-4)  give  the  threshold  condition  as  3.9  ^u)  /u)  ^5.8. 

v  pe  ce 

implies  that 

,  XT  2,1/2 
4ttN  e  \ 

pe  I  ^  eB 


m 


me 


This 


(2-10) 


or, 


N  ^25 
pe 


4iTmc 


(2-11) 


For  B=0.5  gauss.  Equation  2-11  yields  the  critical  plasma  electron  density 


pe 


6  xio5  (cm  3) 


(2-12) 


On  the  other  hand,  Winckler  (Reference  2-18)  suggests  that  the 

2  2 

discharge  may  occur  when  the  Brillouin  limit  (o)^=0.5  co^)  is  reached  where 
refers  to  the  beam  itself.  This  would  imply  that  the  critical  density 
of  beam  electrons,  N  ,  required  for  onset  of  the  discharge  is 


N. 


B" 


be  2 
87imc 


(2-13) 


or,  for  a  field  of  0.5  gauss, 
N. 


be 


=  1.2  x  104  (cm'3) 


(2-14) 


The  optical  consequences  of  a  BPD,  as  measured  in  the  Bernstein 
experiments  (Reference  2-5),  are  that  there  is  at  least  a  10  fold  enhance- 

o 

ment  in  the  intensity  of  3914  A  emission  and  an  increase  in  the  beam  radius 
by  a  factor  of  3.  The  phenomenon  is  also  accompanied  by  strong  radiowave 
emissions.  However,  these  emissions  will  not  be  considered  here  since  no 
instrumentation  to  detect  them  was  present  during  the  EXCEDE  experiments. 
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Relevance  to  a  Nuclear  Environment 


Electron  energy  deposition  in  a  nuclear  environment  occurs  mainly 
through  X-ray  photoelectrons  and  beta  electrons  from  fission  debris.  Are 
either  of  these  processes  likely  to  induce  a  BPD  condition?  Consider  each 
separately : 

X  Ray  Absorption.  In  X-ray  absorption,  the  photoelectrons  are 
emitted  preferentially  at  right  angles  to  the  direction  of  the  incident 
X-ray  photons,  which  move  out  radially  from  the  burst  point.  Under  such 
conditions,  no  electron  beam  is  formed;  consequently,  no  BPD  can  occur. 


Beta-Particle  Streamers.  The  current  density  of  beta  electrons 

2 

can  be  estimated  as  follows.  Let  A  be  the  area  (km  )  over  which  radio¬ 
active  debris  is  uniformly  distributed.  It  can  then  be  shown  (Reference 

2-19)  that  the  number  of  betas,  N0 ,  crossing  unit  area  in  unit  time  is 

p 

given  by 


8.8x  1015  IV 


nb 


A(l+t) 


1.15 


f  ,u  *  -2  -K 

—  (betas  cm  sec  ) 


(2-15) 


where  is  the  fission  yield  (MT) ,  and  t  is  the  time  (sec)  after  burst. 

The  current  density,  JD ,  is  then 

p 


1  o  o 

J6  =  Np/6-25x  10  (amps  cm  )  •  (2-16) 


Furthermore,  if  v^  is  the  beta  particle  velocity,  their  number  density 
is 


nB  =  NVVB  •  (2-17) 

For  a  circular  distribution  of  debris  over  a  radius  R  (km),  the 
value  of  Ng,  per  MT  of  fission  yield,  at  t  =  1  sec,  is 

XT  1.26x10^  .  -2  -1  u  - 1 .  rn  -i  on 

N0(l)  =  - r -  (betas  cm  sec  MT  )  .  (2-18) 

P  2- 
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Equations  2-16,  2-17,  and  2-18  then  yield 


V15  = 


2. Ox  10 

R2 


(amps  cm  2) 


(2-19) 


and,  for  1.5  MeV  beta  particles  with  a  velocity  of  2.6x  10^(1  cm  sec"* , 


nBU)  - 


4 . 9  x  10 
R2 


,  -3. 

(cm  ) 


(2-20) 


Table  2-1  shows  values  of  n^(l)  and  Jg(l)  for  selected  values  of  R. 


Table  2-1.  Beta  particle  concentration  and  current  density 
per  MT  fission  yield  at  1  sec  for  1.5  MeV  betas. 


Debris  Radius 

(km) 

nB0) 

(electrons  cm  ) 

°e  < 1  > 

_  o 

(amps  cm  ) 

10 

4. 9  x  102 

2  x  10"6 

102 

4.9 

2  x  10"8 

103 

4. 9  x  10"2 

2  x  10"10 

104 

4. 9  x  IQ"4 

2  x  )0"12 

For  high-altitude  nuclear  bursts  involving  fission  yields  as  large 

as  1  MT,  the  debris  radius  is  expected  to  exceed  100  km.  Thus,  from  Table 

-3  -8 

2-1,  we  expect  that  in  beta  streamers  nQ  <  5  electrons  cm  '  and  JD  <  2  x  10 
2  p  p 

amps  cm  .  These  values  are  to  be  compared  with  the  corresponding  values, 

suggested  by  laboratory  data,  and  quoted  in  the  previous  subsection,  of 

about  10^  electrons  cm  ^  and  5 x  10  ^  amps  cm  ^  for  the  minimum  beam  electron 

and  current  densities  required  for  onset  of  a  BPD.  It  thus  appears  that  the 
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current  density,  or  the  beam  electron  density  in  beta  streamers  is  much 

★ 

too  low  to  induce  any  BPD  effects. 

We  conclude  that  the  BPD  phenomenon  should  not  be  present  in  a 
nuclear  environment. 

Relevance  to  the  Aurora 


Consider  the  case  of  an  unusually  strong  class  III  aurora  with  a 
primary  electron  flux  at  the  top  of  the  atmosphere  of  the  form 


-E,  y/ 8.5 

$  =  2.2  x  lO^e  C  (electrons  cm  "  sec  eY  )  .  (2-21) 

14  -2  -1 

This  corresponds  to  an  integrated  deposition  rate  of  1.6x  10  eV  cm  sec 

o 

that  results  in  an  auroral  intensity  at  5914  A  of  about  230  kR. 


It  is  easily  shown  that,  for  this  strong  aurora,  the  electron 

-3  -9  -2 

and  current  densities  have  the  values  of  6  cm  and  3  x  10  amps  cm  , 

respectively.  These  numbers  are  clearly  small  compared  with  the  corres¬ 
ponding  values  of  10^  cm  ^  and  5x  10  ^  amps  cm  ^ ,  cited  above,  that  appear 

to  be  required  for  onset  of  a  BPD.  We  thus  conclude  that  the  BPD  phenomenon 

* 

is  not  present  in  an  auroral  environment. 

Relevance  to  EXCEDE  Experiments 


The  average  current  density,  ,  and  beam-electron  density,  , 
for  a  beam  radius  (cm),  are  given,  respectively,  by 

I, 


Jb  = 


rV 


(amps  cm  ^) 


(2-22) 


*  On  the  other  hand,  the  minimum  necessary  plasma  density  requirement, 
stated  by  Equation  2-12,  is  easily  satisfied  in  a  beta  patch  or  in  a 
strong  auroral  arc.  To  our  knowledge,  no  BPD  effects  have  been  observed 
in  the  aurora. 
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and 


1  o 

6.25  x  10  J 

- - - -  (cm"3)  ,  (2-23) 

where  1^  is  the  beam  current  (amps),  and  v  is  the  velocity  ( cm  sec"1) 
of  the  beam  electrons. 

In  the  EXCEDE : SPECTRAL  experiment,  the  nominal  electron  energy 
9  -1 

was  3  keV  (v=2. 29x10  cm  sec  ),  and  the  current  was  7  amps  per  electron 

gun.  Table  2-2  shows  values  in  this  case  for  J,  and  N,  ,  inferred  from 

b  be 

Equations  2-22  and  2-23,  for  selected  values  of  the  beam  radius. 


Table  2-2.  Average  beam-electron  and  current  densities  in 
the  EXCEDErSPECTRAL  experiment. 


Beam  Radius 

(m) 

N. 

be 

-3 

(electrons  cm  ) 

Jb 

_  o 

(amps  cm  ) 

1 

6.1  x  105 

2. 2  x  10-4 

2 

1 . 5  x  1 05 

5. 6  x  10"5 

3 

6.8  x  104 

2. 5  x  10"5 

4 

o 

X 

00 

CO 

1 .4x  10"5 

5 

o 

X 

(XI 

8.9  x  io"6 

The  effective  beam  radius,  inferred  from  the  onboard  camera  data 

(Reference  2-9),  was  about  2  m.  For  this  radius,  we  see  from  Table  2-2 

that  the  values  for  N^e  and  exceed  the  corresponding  critical  values 

4  "  / 

(of  about  10  and  5x10  ° ,  respectively)  for  BPD  onset  by  at  least  an  order 
of  magnitude.  Furthermore,  to  the  extent  that  it  is  relevant,  the  plasma 
density  at  123-km  altitude  in  EXCEDErSPECTRAL  builds  up  very  rapidly  to 
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values  greatly  exceeding  the  critical  value  of  6  x  10  given  by  Equation 
2-12.  Consequently,  all  of  the  criteria  listed  above  for  onset  of  a  BPD 
seem  to  be  well  satisfied  by  EXCEDE  electron  beams. 

In  addition,  the  optical  consequences  of  a  BPD,  as  inferred  from 
the  Bernstein  experiments,  also  seem  to  be  satisfied,  at  least  for  the 
EXCEDE : SPECTRAL  image  we  have  studied.  In  particular,  the  laboratory - 

o 

observed  enhancement  in  the  3914  A  intensity  by  at  least  an  order  of  magni¬ 
tude  is  certainly  well  fulfilled  by  the  EXCEDE  results.  Furthermore,  the 
laboratory-observed  increase  in  the  beam  radius  by  a  factor  of  3,  after  BPD 
onset,  also  appears  to  be  duplicated  in  EXCEDE.  This  can  be  seen  by  com¬ 
paring  the  dashed  and  upper  solid  curves  in  Figure  2-2.  The  upper  solid 
curve  is  the  calculated  radial  intensity  distribution,  increased  by  a  factor 
of  10,  and  shifted  to  bring  its  peak  into  line  with  that  of  the  dashed  curve, 
which  represents  the  measured  radial  distribution  4  meters  in  front  of  the 
electron  gun.  The  full  width  at  half  maximum  of  the  dashed  curve  is  almost 
exactly  a  factor  of  3  greater  than  that  of  the  solid  curve. 

The  conclusion  then  is,  that  the  laboratory  conditions  for  onset 
of  a  BPD,  as  well  as  certain  optical  properties  of  the  resulting  emission, 
appear  to  have  been  fulfilled  in  the  EXCEDE : SPECTRAL  experiment. 
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SECTION  3 


UTILITY  OF  EXCEDE  FOR  PLASMA 
RADIATION  VERIFICATION 


PRELIMINARY 

It  has  been  established  (Reference  1-1)  that  plasma  radiation 
from  the  fireball  plume  of  a  high-altitude  nuclear  detonation  is  the  single 
greatest  potential  source  of  interference  with  LWIR  sensor  systems.  This 
radiation  includes  the  continuum  contribution  resulting  from  free-free  and 
free-bound  collisions  between  electrons  and  ions,  as  well  as  line  radiation 
from  bound-bound  transitions  that  occur  following  radiative-col lisional  re¬ 
combination  in  an  0+/N+  plasma.  Although  the  formula  used  to  describe  the 
free-free+free-bound  contribution  is  believed  to  be  on  solid  theoretical 
ground,  the  models  that  describe  the  bound-bound  radiation  have  not  been 
adequately  verified  by  the  limited  experiments  (References  3-1,  3-2) 
performed . 


Spectroscopic  observations  of  the  late-time  Kingfish  and  Checkmate 
fireballs  indicate  that,  during  radiative-col lisional  recombination,  many 
of  the  excited  states  of  the  neutral  atoms  from  which  photon  emission  arises 
are  not  in  local  thermodynamic  equilibrium  (LTE)  with  the  ground  state  of 
the  atom  or  of  the  parent  ion  (Reference  3-3) .  This  means  that  the  state 
populations  cannot  be  determined  simply  from  the  Boltzmann  Equation,  or 
the  Saha  Equation,  but  rather  from  detailed  non  LTE  calculations  involving 
the  rates  for  electron-ion  recombination  and  for  collisional  and  radiative 
cascade  to  the  ground  state  of  the  atom.  Many  uncertainties  are  involved 
in  these  non  LTE  calculations;  consequently,  experimental  verification  of 
the  results  is  highly  desirable. 


41 


The  experimental  requirements  for  verification  of  the  bound-bound 
radiation  models,  under  non  LTE  conditions,  have  been  reviewed  by  Sappenfield 
(Reference  3-3).  The  principal  requirement  is  that  the  plasma  density  be 
reasonably  high  while,  at  the  same  time,  the  electron  temperature  be  low. 

If  both  quantities  are  large,  the  excitation  rate  by  collisions  may  be 
sufficiently  great  to  assure  LTE  conditions  which,  for  the  experiment,  we 
wish  to  avoid. 

The  question  arises  as  to  whether  the  requirements  for  verifica¬ 
tion  might  be  met  by  an  EXCEDE-type  experiment.  To  help  answer  this  question 
we  made  computer  calculations  using  the  ARCTIC  code.  The  results  of  that 
study  are  reported  in  this  section. 

As  explained  in  Section  2,  the  ARCTIC  code  was  originally  designed 
for  use  under  auroral  bombardment  conditions  that  are  about  10^  times  less 
intense  than  those  encountered  in  EXCEDE  experiments.  As  a  consequence, 
considerable  time  and  effort  was  required  to  modify  the  code  so  that  it 
would  be  applicable  under  these  extreme  conditions,  particularly  with  respect 
to  a  determination  of  the  electron  temperature.  A  brief  outline  of  the 
modifications  that  were  made  in  ARCTIC  are  given  in  Section  2,  and  a  more 
detailed  description  can  be  found  in  Appendix  A. 

ARCTIC  CODE  CALCULATIONS 

The  calculations  have  been  done  for  the  case  of  EXCEDE : SPECTRAL, 

at  altitudes  of  90  and  123  km.  In  both  cases  we  have  assumed  the  maximum 

energy  deposition  rate,  which  occurs  near  the  electron  gun,  as  determined 

by  our  earlier  code  calculations.  This  rate,  e,  in  normalized  units  for  a 

3  -1  -1  -I 

single  electron,  is  1.97*10"  eV  gm  sec  electron  (Reference  2-1, 

•  —  ^  ^ 

Figure  3-3).  The  deposition  rate,  in  units  of  eV  cm  '  sec  ,  is  then 
given  by 

<?  =  6. 25  x  1018  Ipe  ,  (3-1) 
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where  I  is  the  beam  current  (amps),  and  p  is  the  atmospheric  density 
-  3 

(gm  cm  ).  For  a  7  amp  beam,  the  deposition  rate  is  then 

&  =  8.62  x  lo22p  .  (3-2) 

Species  Concentrations  and  the  Electron  Temperature 

-11  -3 

123-km  Altitude.  From  Equation  3-2,  with  p=  1.72*  10  gm  cm  , 

.  12  -3  -1 

the  energy  deposition  rate  input  to  the  code  is  1.48x  10  eV  cm  sec 

_7 

The  code  was  run  for  times  from  10  sec  to  10  sec.  Some  of  the  results 
are  shown  in  Figure  3-1.  Included  are  values  for  selected  species  concen¬ 
trations  as  well  as  the  electron  temperature. 

Some  points  to  notice  from  Figure  3-1  are  the  following.  The 

electron  density  builds  up  rather  uniformly  until  about  10  ^  sec  when  it 

begins  to  level  off  because  it  is  approaching  its  quasi-equilibrium  value 
9  -3 

(^1.5*  10  cm  '  )  ,  controlled  by  dissociative  recombination  with  the 
molecular  ions,  which  tend  to  dominate  the  ion  concentration  until  about 
1  sec.  However,  dissociation  of  and  0^  by  the  continual  electron  bombard¬ 
ment  results  in  a  buildup  in  the  N,  N+,  0,  and  0+  concentrations  to  the 

point  where  the  ionization  is  dominated  by  the  species  N  and  0  after  about 

1.5  sec.  Recombination  is  then  controlled  by  the  slow  radiative-col 1 isonal 
process  so  that  the  electron  density,  along  with  the  N+  and  the  0+  concentra¬ 
tion,  rises  rapidly. 

The  electron  temperature,  ,  first  rises  uniformly  because  the 

heating  rate  by  the  bombarding  flux  dominates  the  cooling  rate  from  all 

-4 

processes  until  about  2.5 x  10  sec  when  T^  attains  its  maximum  value  of 
about  45,000  °K  and  the  two  rates  are  then  equal.  Subsequently,  until 

about  10  1  sec,  cooling  by  electronic  and  vibrational  excitation  of 

dominates  the  heating,  and  the  temperature  drops  to  about  6000  °K,  a  value 
too  low  to  maintain  significant  excitation  of  the  At  this  point,  heating 

by  the  bombarding  flux  and  by  deexcitation  of  the  metastable  states  of  0+, 
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Figure  3-1.  Calculated  concentrations  of  selected  species,  and  the  electron  temperature,  at 
123-km  altitude  under  EXCEDE  bombardment  conditions.  (Energy  deposition  rate  = 


N+,  N  and  0  dominates  the  cooling,  and  then  rises  until  about  2  sec 

when  the  heating  and  cooling  rates  are  equal.  However,  the  rapid  rise  in  the 
atomic  ion  density  after  1.5  sec  produces  a  precipitously  large  rate  of 
electron  cooling  by  elastic  collisions  with  the  ions,  and  T  drops  rapidly 

A  ^ 

from  a  value  of  about  104  °K  at  3.6  sec  to  about  400  °K  at  4  sec.  Accom¬ 
panying  this  drop  in  T  is  a  corresponding  increase  in  the  value  of  the 

G  - 14  3  - 1 

radi  ative-col  li  sional  recombination  coefficient  from  6.2x10  cm  sec 

at  3.6  sec  to  1.3*  10  ^  cm‘^  sec  ^  at  4  sec.  This  produces  a  sudden  drop 

in  the  electron  and  ion  densities,  and  thus  in  the  electron  cooling  rate  by 

elastic  collisions  with  the  ions,  so  that  the  electron  temperature  recovers 

somewhat  and  by  5  sec  attains  an  equilibrium  value  of  about  500  °K,  while 

the  electron  density  is  close  to  its  new  equilibrium  value  of  about 

1.6x  10^°  cm  By  10  sec,  about  75  percent  of  the  ions  are  N+;  the 

remainder  are  0+. 

-9  -3 

90-km  Altitude.  From  Equation  3-2,  with  p  =  3.4x  10  gm  cm  , 

— —  ■  _  3  _  i 

the  energy  deposition  rate  input  to  the  code  is  2.93x  10  eV  cm  sec 

.7 

Results  were  obtained  for  times  from  10  sec  to  5.8  sec.  Some  of  these, 
including  values  for  selected  species  concentrations  and  the  electron 
temperature,  are  shown  in  Figure  3-2. 

The  commentary  on  Figure  3-2  generally  follows  that  for  Figure  3-1. 

The  electron  density  bui  lds  up  rather  uniformly  until  about  2  x  10  sec  when 

it  starts  to  level  off  and  then  attains  its  quasi-equilibrium  value  of  about 
10  -3 

1.2  x  10  cm  which  persists  until  about  1.5  sec.  During  this  time  the 
molecular  ions  N0+,  0n+,  and  NO+  are  dominant  so  that  the  quasi-equilibrium 
electron  density  is  controlled  by  dissociative  recombination.  However, 
dissociation  of  and  0?  by  the  bombarding  flux  results  in  a  buildup  in 

■f  ^ 

the  N,  N  ,  0,  and  0  concentrations  so  that  after  1.5  sec  the  ionization  is 

dominated  by  N+  and  0+  ions.  Radiative-col lisional  recombination  then  takes 

over  from  dissociative  recombination,  and  the  electron  density,  along  with 

the  N+  and  the  0+  concentration,  rises  rapidly  and  by  5.8  sec  the  electron 

11  -3 

density  attains  a  value  of  about  3  x 10  cm  .  The  code  was  not  pushed 
beyond  5.8  sec  because  of  excessively  long  running  times. 
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Time  (sec) 

Figure  3-2.  Calculated  concentrations  of  selected  species,  and  the  electron  temperature,  at 
90-km  altitude  under  EXCEDE  bombardment  conditions.  (Energy  deposition  rate  = 
2.93  x  10*4  ey  cm-3  Sec-1 ) 


Until  about  1  sec,  the  electron  temperature  behavior  is  approxi¬ 
mately  similar  to  that  described  above  in  the  123-km  case  for  times  to  2  sec. 

From  a  value  of  about  6000°K  at  1  sec,  T  then  gradually  subsides,  mainly 

1  1  C  3 

by  excitation  of  0(  D)  and  0(  S)  and  fine  structure  excitation  of  0(  P)  , 

until  about  2.5  sec  when  the  equilibrium  (or  quasi  equilibrium)  temperature 

of  about  4500°K  is  attained.  The  rapid  drop  in  seen  above  in  the  123-km 

case  has  not  occurred  in  this  case,  at  least  by  5.8  sec.  This  may  be  because 

the  late-time  cooling  mechanisms  are  different  in  the  two  cases.  At  123  km 

the  electron  cooling  is  mainly  by  elastic  collisions  with  the  ions;  at  90  km, 

it  is  mainly  by  excitation  of  atomic  oxygen. 

Plasma  Radiation 


The  volume  emission  rate  from  free-free  (ff)  +  free-bound  (fb) 
collisions  of  electrons  with  ions  can  be  written  as  the  frequency  (v)-inde- 
pendent  expression  (Reference  3-4) 


J  (ff+fb)  = 

V  3/6 


8, 1/2  /X 


6 

e 

3/2  3 


N  N 
+  e 


m 


(kTe) 


172 


,  -3  - 1  -1  u  -K 

(ergs  cm  sec  ster  Hz  ) 

(3-3) 


Here,  e  and  m  are  the  electron  charge  and  mass,  respectively,  and 

N  are  the  electron  and  positive  ion  concentrations,  respectively,  is 

the  electron  temperature  (°K),  k  is  Boltzmann's  constant,  and  g/g+  is 
the  ratio  of  statistical  weights  of  the  ground  states  of  the  neutral  and 
ion  species.  For  hydrogen-like  ions,  which  we  assume  here,  g/g+=2. 
Alternatively,  with  a  change  of  units,  Equation  3-3  becomes 

-3  -1  -1 

(watts  cm  ster  Hz  )  (3-4) 

-3  -1  -1 

(watts  cm  ster  pm  )  (3-5) 

where  A  is  the  wavelength  in  micrometers. 


J  (ff+fb)  =  5 . 39  x  10"46N  N  T  '1/“ 
v v  +  e  e 


or 


J,  (ff+fb)  =  1.62  x  10 
P 


N  N  T 
■31  +  e  e 


1/2 


X, 
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For  the  volume  emission  rate  from  bound-bound  (bb)  transitions 
we  have  used  the  approximate  expression  (Reference  3-5) 


J.  (bb)  =  2.65  x  10 

A 


-20 


ctN  N 
+  e 


T  °-2'45  ‘O^IOV1-97 
e  U 


(watts  cm  ^  ster  *  ym  ^)  (3-6) 


where  a,  the  radiative-col lisional  recombination  coefficient,  is  deter¬ 
mined  by 


a  ^  5. 75  x  10 


•13N 

e  \ 11606/ 


[-0.0501(log10Ne)  +0. 4641  1 °g10Ne“ 1 ‘ 913255] 


(3-7) 


and  T  is  the  neutral  termperature  (°K).  Although  Equation  3-6  is  a  very 
low  resolution  approximation  that  does  not  describe  the  spikes  associated 
with  individual  spectral  lines,  we  believe  that  it  is  adequate  for  ev¬ 
aluating  EXCEDE's  utility  as  a  plasma  radiation  experiment. 


Equations  3-5  through  3-7  were  used,  together  with  values  of 
N+ ( [0+] + [N+] ) ,  N^,  and  from  Figures  3-1  and  3-2,  to  obtain  plots  of 

the  quantity  [J^  (ff+fb) +J^  (bb)  ] .  In  these  calculations,  neutral 

temperatures  of  175.0  and  367.3  °K  were  used  for  altitudes  of  90  and  123  km, 
respectively.  The  results  are  shown  in  Figure  3-3. 


For  the  123-km  case,  the  calculations  show  that  the  free-free+ 
free-bound  contribution  exceeds  the  bound-bound  contribution,  for  times  to 
3.8  sec,  by  factors  that  vary  between  about  1.6  and  5.6.  However,  for  times 
after  3.8  sec,  the  bound-bound  contribution  is  dominant  by  at  least  2  orders- 
of-magnitude.  This  is  because  of  the  rapid  drop  in  the  electron  temperature 
at  this  time  that  results  in  a  large  increase  in  the  radiative-collisional 
recombination  coefficient,  as  described  above. 


For  the  90-kmcase,  the  bound-bound  contribution  exceeds  the  free¬ 
ly 

free+free-bound  contribution  at  all  times  from  about  10  sec  to  5.8  sec 
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Figure  3-3.  Calculated  volume  emission  rate  from  plasma  radiation 
(ff+fb+bb)  at  90-  and  123-km  altitude  in  an  EXCEDE- 
type  experiment. 
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by  factors  that  vary  between  about  1.2  and  4.  At  about  1.5  sec,  the  rapid 
rise  in  the  electron  density  produces,  by  about  5.7  sec,  a  4  order-of- 
magnitude  increase  in  the  plasma  emission  rate.  However,  because  is 

still  fairly  large,  the  bound-bound  contribution  at  that  time  is  only  a  fac¬ 
tor  of  4  larger  than  the  free-free  +  free-bound  contribution. 

APPLICATION  TO  EXCEDE 

As  far  as  EXCEDE  is  concerned,  the  main  point  to  note  about  both 
curves  in  Figure  3-3  is  that  the  rapid  increase  in  the  emission  rate  does 
not  occur  until  after  one  or  more  seconds  of  bombardment  time.  This  is  long 
compared  with  typical  bombardment  times  in  EXCEDE  experiments  that  are 
usually  limited  to  less  than  about  0.1  sec  because  of  rocket  motion.  For 
bombardment  times  this  short,  it  is  clear  from  our  calculations  that  not 
only  is  the  magnitude  of  the  plasma  emission  small,  but  also  the  bound- 
bound  contribution,  that  we  are  especially  interested  in  measuring,  is 
either  less  than  the  free-free+free-bound  continuum  (123-km  case)  or  is  not 
significantly  larger  than  it  (90-km  case). 

To  achieve  bombardment  times  of  a  few  seconds  duration  in  an 
EXCEDE  experiment  would  require  flying  the  rocket  parallel  to  the  earth's 
magnetic  field  for  a  similar  length  of  time.  Although  we  are  not  familiar 
with  the  logistics  involved  in  operating  the  rocket  in  this  manner,  we 
suspect  that  they  are  sufficiently  difficult  to  render  that  mode  of  opera¬ 
tion  infeasible. 

At  this  point  we  should  emphasize  that  the  foregoing  results  and 
considerations  are  based  on  our  code  calculations  that  assume  the  validity 
of  the  independent  particle  model  in  which  the  excitation  and  ionization  is 
produced  directly  by  the  beam  electrons.  However,  as  discussed  in  Section  2, 
this  model  fails  to  account  for  the  magnitude  and  dimensions  of  the  intense 
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optical  emission  observed  in  the  vicinity  of  the  EXCEDE : SPECTRAL  guns.  It 
has  been  suggested  that  a  BPD  may,  in  fact,  be  responsible  because  the  con¬ 
ditions  for  its  ignition  appear  to  have  been  satisfied  in  the  EXCEDE  experi¬ 
ments.  If  such  is  indeed  the  case,  then  the  behavior  of  the  plasma  radiation, 
with  respect  to  both  its  build-up  time  and  its  magnitude,  may  be  quite  dif¬ 
ferent  from  that  shown  in  Figure  3-3.  For  example,  it  has  been  observed 
that  BPD's  are  ignited  within  a  few  microseconds  after  beam  turn  on  and 
that  plasma  electrons,  not  beam  electrons,  are  then  the  dominant  ionization 
source  (Reference  2-3).  Under  these  conditions,  one  would  expect  a  much 
earlier  buildup  in  the  plasma  density  and,  possibly,  in  the  intensity  of 
plasma  radiation.  The  argument  cited  earlier  against  EXCEDE  as  a  useful 
tool  for  performing  a  plasma  radiation  experiment,  based  on  bombardment¬ 
time  considerations,  would  then  not  be  valid. 

On  the  other  hand,  if  BPD  conditions  in  EXCEDE  should  "save  the 
day"  in  terms  of  generating  large  bound-bound  emission  rates  at  early  times, 
there  would  still  be  the  problem  of  making  accurate  measurements  of  the 
electron  temperature  and  atomic  ion  concentrations  in  the  radiating  region. 

It  seems  to  us  that  an  experiment  of  this  nature  could  be  done  more  readily, 
and  probably  less  expensively,  in  the  laboratory. 

COMPETING  MOLECULAR  RADIATIONS 

We  have  compared  the  magnitude  of  the  calculated  plasma  radiation, 
integrated  over  the  8-  to  12-ym  and  2-  to  5-ym  bands,  with  competing  molecular 
emission  over  the  same  wavelength  regions.  The  results  are  shown  below. 

For  the  high  deposition  rates  inherent  in  EXCEDE,  we  find  that  the 
IR  molecular  emission  is  completely  dominated  by  N0  fluorescence  bands. 

*  However,  if  the  effective  electron  temperature  is  very  high  under  BPD  con¬ 
ditions,  because  of  internal  electric-field  heating,  then  the  bound-bound 
plasma  emission  rate  may,  in  fact,  not  be  large. 
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An  estimate  of  this  fluorescence  was  made  using  some  (unpublished)  results 
obtained  from  the  FLUOR  code  (RDA)  and  DANTE  code  (TEMPO)  (Reference  3-6) 
in  which  we  increased  the  Wu-Benesch  contribution  by  a  factor  of  20  to  con¬ 
form  with  optical  data  from  the  EXCEDE : SPECTRAL  experiment  (Reference  3-7). 
Figures  3-4  and  3-5  show  the  resulting  spectral  efficiency,  n~,  in  units 
of  ergs  per  cm  per  100  photons  at  3914  A,  in  the  region  from  200  to  5000 
cm  \  for  altitudes  of  92.5  and  120  km,  respectively.  The  fraction,  f^, 
of  the  deposited  energy  that  is  radiated  over  the  wavenumber  interval 
V1  to  V2  Can  °^taine(i  from  the  relation 


1 


a2 


8.9x10  j  o  ~  dv 


(3-8) 


From  the  previous  subsection,  the  energy  deposition  rates  assumed 

for  altitudes  of  90  and  123  km  are  2.93x10^  and  1.48x10^  eV  cm  ^  sec  \ 

respectively.  These  values,  together  with  the  fluorescence  efficiency 

curves  in  Figures  3-4  and  3-5,  lead  to  the  volume  emission  rates  for 

molecular  radiation  shown  in  Table  3-1.  For  both  wavelength  regions,  the 

3  3 

principal  contributors  are  the  Wu-Benesch  (W  A^-B  11^)  and  First  Positive 

(B^IT  -A^Z+)  band  systems  of  N_  . 
g  ir  2 

Table  3-1.  Comparison  between  plasma  and  molecular  radiation. 


A1 ti tude 

Wavelength  Band 

Volume  Emission  Rate 

-  3  - 1 

(watts  cm  ster  ) 

(km) 

(ym) 

Molecular  Radiation 

Max.  Plasma  Radiation 

90 

2-5 

1.1  x  10"8 

O 

1 

o 

X 

o 

CO 

8-14 

4. 8  x  10"11 

5.4  x  10'11 

123 

2-5 

7.8  x  10"11 

1 . 2  x  10"10 

8-14 

-1  3 

3. 8  x  10 

2.1  x  10"11 
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Figure  3-4b.  N?  fluorescence  efficiency  at  92.5  km  (T  =  192°K). 
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Figure  3-5a.  N?  fluorescence  efficiency  at  120--km  altitude  (T  =  355°K). 
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Figure  3-5b.  N9  fluorescence  efficiency  at  120-km  altitude  (T  =  355°K). 


To  compare  the  molecular  radiation  with  the  corresponding  plasma 
radiation,  we  note  that  since  the  latter  varies  approximately  inversely  as 
\j  (Equations  3-5  and  3-6),  the  plasma  emission,  integrated  between  wave¬ 
lengths  \  and  X2,  is 

X  X2 

f  dX  /  i  i  \ 

dX^  =  K  /  — ~  j —  -  —  J  K  (watts  cm  ster  )  (3-9) 


where  K  is  plotted  in  Figure  3-3.  The  maximum  values  of  K  shown  in  this 

-9  -10  -3  -1  -1  2 

figure  are  approximately  1x10  and  5x10  watts  cm  ster  pm  (dm) 

for  the  90-km  and  123-km  cases,  respectively.  Use  of  these  values  in  Equa¬ 
tion  3-9  then  leads  to  the  plasma  emission  rates  shown  in  Table  3-1. 

Points  to  note  about  the  numbers  shown  in  Table  3-1  are  the  fol¬ 
lowing.  In  the  90-km  case,  it  appears  that  the  molecular  emission  dominates 
the  plasma  radiation  over  the  2-5  pm  band,  and  that  the  two  are  nearly  equal 
over  the  8-14  pm  band.  In  the  123-km  case,  the  maximum  plasma  emission  rate 
(almost  exclusively  from  bound-bound  radiation)  exceeds  the  molecular  emis¬ 
sion  by  nearly  2  orders-of-magnitude  in  the  LWIR  band  (8-12  pm),  but  by  a 
much  smaller  factor  in  the  SWIR  band  (2-5  pm).  Thus,  in  the  context  of  a 
plasma  radiation  experiment,  the  signal  to  noise  ratio  in  the  LWIR  is  sig¬ 
nificantly  greater  than  it  is  at  the  shorter  wavelengths. 


In  conclusion,  it  appears  that  molecular  radiation,  induced  by 
electron-beam  excitation  of  air,  should  not  seriously  interfere  with  observa¬ 
tions  of  plasma  radiation  in  the  LWIR  region  at  altitudes  as  high  as  123  km 
once  the  plasma  density  has  been  built  up  sufficiently  for  electron-ion  col¬ 
lisions  to  cool  the  plasma.  However,  this  conclusion  does  not  necessarily 
hold  up  under  beam  plasma  discharge  conditions,  particularly  if  the  internal 
electric  fields  should  maintain  a  high  value  for  the  electron  temperature 
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that  would  then  serve  to  depress  the  radiative-col  1 isional 
efficient  and,  hence,  the  bound-bound  plasma  emission  rate 
under  such  conditions,  collisional  excitation  of  radiation 
electrons  might  well  mask  the  recombination  radiation. 


recombination  co- 
Furthermore , 
by  the  hot  plasma 
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SECTION  4 

SOME  LWIR  CONSIDERATIONS 


This  section  describes,  very  briefly,  three  LWIR- related  topics 
that  we  considered  as  part  of  the  work  statement  under  the  present  contract. 
They  are:  (1)  the  possible  excitation  of  LWIR  emission  from  N\0  by  solar 
pumping,  (2)  photon  excitation  of  the  LWIR  bands  of  uranium-  and  fission- 
fragment  oxide  ions,  and  (3)  preliminary  assessment  of  how  the  relative 
importance  of  LWIR  emission  mechanisms,  inferred  from  single-burst  calcu¬ 
lations,  may  be  modified  in  a  multiburst  environment. 

The  need  to  consider  Topics  1  and  2  was  discussed  in  Reference 
1-1  which  treats  the  relative  importance  of  various  LWIR  emission  mechan¬ 
isms  in  a  single-burst  environment.  Our  present  work,  as  related  to  these 
two  topics,  has  consisted  primarily  of  literature  surveys.  Item  3  above 
is  a  natural  outgrowth  of  the  work  reported  in  Reference  1-1.  As  yet,  no 
in-depth  study  has  been  possible  for  a  multiburst  environment  due  to  the 
very  limited  amount  of  calculational  data  available.  This  situation  should 
be  rectified  in  the  near  future  when  the  optical/lR  capability  of  the 
SCENARIO  code  is  fully  operative. 

SOLAR  PUMPING  OF  LWIR  BANDS  OF  N20 

The  mechanism  involved  here  is  one  in  which  a  molecule  is  first 
electronically  excited  by  UV-  or  visible-light  absorption  in  its  ground 
state.  Subsequently,  following  radiative  cascade,  it  may  be  left  in  a 
vibrationally  excited  state  of  the  ground  electronic  state.  If  the 
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molecule  is  IR  active,  further  radiative  cascade  will  result  in  the  emission 
of  its  characteristic  IR  bands. 


It  can  be  shown  (Reference  4-1,  p.  382)  that  the  emission  rate, 
cp$>  of  IR  photons  per  molecule  due  to  solar  pumping  is  determined  from  the 
equation 


Tre  _  - 
me  ^  e 


(4-1) 


Here,  p  is  the  number  of  IR  photons  emitted  per  visible/UV  photon  absorbed, 
is  the  solar  photon  flux  (photons  cm  ^  sec  Hz  ),  and  f^_  is  the 
oscillator  strength  for  the  electronic  transition  involved.  Clearly,  two 
requirements  for  strong  IR  emission  by  this  mechanism  are:  first,  that 
there  be  strong  absorption  of  visible  or  UV  photons  by  the  ground  state  into 
an  upper  state  that  can  radiate  back  to  the  ground  (electronic)  state  again 
(i.e.  large  f^) ;  second,  that  the  solar  flux,  ,jF^,  be  relatively  large. 


The  first  requirement  precludes  situations  in  which  the  absorption 
(strong  though  it  may  be)  is  to  an  upper  state  that  spontaneously  dissociates, 
or  to  one  that  overlaps  a  continuum  so  that  predissociation  (or  preioniza¬ 
tion)  may  occur.  The  case  of  predissociation  may  result  in  diffuse  visible/ 
UV  band  radiation,  but  it  is  not  likely  to  produce  strong  IR  emission.  The 

second  requirement  favors  absorption  by  visible- light  photons,  rather  than 

3 

UV  photons,  because  the  solar  flux  varies  approximately  as  A  .  Thus,  for 

o 

example,  other  things  being  equal,  absorption  of  5500  A  photons  will  be 
about  60  times  more  effective  in  producing  IR  emission  than  will  the 

o 

absorption  of  1400  A  photons. 


We  have  surveyed  the  literature  to  determine  the  availability  of 
suitable  states  of  No0  that  may  permit  an  effective  solar-pumping  mechanism 
to  operate.  The  results  are  negative.  In  the  first  place,  No0  does  not 
absorb  in  the  visible  spectral  region  (Reference  4-2).  As  for  the  UV  region, 
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it  is  dominated  by  continuous  absorption  spectra,  superimposed  in  certain 
spectral  intervals  by  diffuse  bands  (Reference  4-3).  Thus,  qualitatively, 
the  findings  violate  both  criteria  stated  above  for  an  effective  solar¬ 
pumping  mechanism. 

Numerically,  one  may  attempt  to  place  an  upper  limit  on  the  IR- 
photon  emission  rate  by  considering  the  measured  total  absorption  cross 

o 

sections  for  NO.  For  example,  in  the  spectral  region  from  about  1400  A 

o  4 

to  1600  A  the  absorption  spectrum  consists  of  an  apparent  continuum  on 

which  are  superimposed  about  20  diffuse  bands.  A  cross  section  peak  of 
- 18  2  ° 

8x10  cm  occurs  at  about  1450  A  (Reference  4-4).  An  integration  of 

-2 

the  published  cross  section  then  leads  to  a  value  for  fe  of  about  2  x  10 

o 

Also,  since  <^,^10  at  1450  A,  Equation  4-1  yields,  for  /rj,  a  value  of 
about  5x10  ^ . 

This  number  for  A]  is  comparable  to  the  corresponding  numbers 
inferred  for  LIVIR  emission  from  N^O  excited  by  earthshine  scatter  (Reference 
1-1,  p,  84).  However,  since  a  major  part  of  the  UV  absorption  cross  section 

appears  to  be  attributable  to  the  continuum,  we  must  conclude  that  the  actual 

•  -  3 

value  of  /p  for  N^O  is  much  smaller  than  5  x  10  .  That  is,  we  conclude 

that  for  N^O  the  solar-pumping  mechanism  is  much  less  effective  in  exciting 

the  LWIR  bands  of  No0  than  is  the  earthshine-scatter  mechanism. 

OSCILLATOR  STRENGTHS  FOR  URANIUM- 
AND  FISSION-FRAGMENT  OXIDE  IONS 

In  Reference  1-1  we  concluded  that,  next  to  plasma  radiation  from 
plume  regions  of  high-altitude  bursts,  the  most  important  source  of  molecular 
radiation  in  a  nuclear  environment  arises  from  uranium-  and  certain  fission- 
fragment  oxide  ions.  The  uranium  species  involved  are  principally  lJ0o+  and 
U0+.  Our  estimate  of  their  importance  was  based  on  assumed  values  for 
excitation  of  their  IR  bands  by  solar  and  earthshine  radiation  equal  to 
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those  for  A£0.  This  excitation  depends  on  the  magnitude  of  the  electronic 
(for  solar  pumping)  and  vibrational  (for  earthshine  scatter)  oscillator 
strengths,  or  f  numbers.  The  need  to  know  the  f  numbers  for  the  above- 
mentioned  ions  is  clearly  apparent  in  view  of  the  potential  these  ions  have 
for  creating  a  structured  LWIR  nuclear  environment  that  may  degrade  the 
performance  of  LWIR  sensor  systems. 

With  the  hope  of  finding  some  experimental  evidence  relating  to 
these  f  numbers,  we  conducted  a  search  of  the  available  literature.  Not 
too  surprisingly,  no  hard  data  were  found. 

Since  that  time,  theoretical  estimates  have  been  reported  (Refer¬ 
ences  4-5,  4-6)  for  the  vibrational  f  number  of  UO  and  U0+.  The  upper- 
limit  estimate  of  5  x  10  ^  for  the  fundamental  bands  of  these  species 
(A  ^12  pm)  leads  to  a  value  of  about  0.2  photons  sec  ^  molecule  \  similar 
to  the  number  assumed  in  Reference  1-1.  Estimates  have  also  been  made  of 
the  solar  pumping  rate,  (p  ,  (see  Equation  4-1)  in  which  values  for  r\ 
between  1  and  10  and  values  for  f  between  0.1  and  1  were  assumed  (Refer- 

ence  4-7).  These  lead  to  numbers  for  cp  ranging  from  about  2  to  200  IR 
- 1  - 1  s 

photons  sec  molecule  ;  numbers  that  are  from  1  to  3  orders-of-magnitude 
larger  than  the  value  adopted  in  Reference  1-1.  However,  these  values  are 
purely  speculative  at  the  present  time,  being  based  on  arguments  related  to 
the  anticipated  density  of  accessible  electronic  states.  Firm  numbers  for 
the  oscillator  strengths,  both  vibrational  and  electronic,  must  await  the 
results  of  successful  laboratory  experiments;  hopefully,  those  currently 
being  conducted  at  the  Los  Alamos  National  Laboratory  (Reference  4-8) . 


*  The  calculated  values  for  the  10.4-pm  band  of  A£0,  reported  in 
Reference  1-1,  are  nearly  the  same  for  both  solar  pumping  and 
for  earthshine  scatter  —  about  0.2  photons  sec_l  molecule"^. 
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MULTIBURST  ENVIRONMENT 


As  mentioned  earlier  in  this  subsection,  the  conclusions  reached 
in  Reference  1-1  regarding  the  relative  importance,  or  unimportance,  of 
various  LWIR  emission  processes  in  a  high-altitude  nuclear  environment  were 
based  on  calculations,  and  other  considerations,  for  a  single-burst  situa¬ 
tion.  However,  since  a  potential  nuclear  conflict  is  almost  certain  to 
involve  multiple  bursts,  it  behooves  us  to  consider  how  the  foregoing  con¬ 
clusions  might  be  modified  in  a  multiburst  situation. 

One  of  the  conclusions  reached,  for  example,  is  that  LWIR  emission 
from  0^,  and  C0^  is  not  significantly  enhanced  above  about  100-km  altitude 
by  a  single  megaton-class  high-altitude  burst.  Is  this  conclusion  still 
valid  in  a  multiburst  environment?  To  provide  an  adequate  answer  to  this 
question,  and  a  confirmation  of  the  other  conclusions,  we  need  to  perform 
multiburst  calculations  using  a  version  of  the  SCENARIO  code  that  incor¬ 
porates  an  IR  capability.  Unfortunately,  this  version  of  SCENARIO,  which 
is  currently  being  developed  at  MRC,  has  not  been  available  for  our  use 
under  the  present  effort. 

We  have  made  a  preliminary  investigation,  however,  using  the 
only  IR-multiburst  results  available  to  us.  These  results  were  based  on 
calculations  performed  at  MRC  a  number  of  years  ago,  using  the  IRCHEM  code, 
for  a  case  in  which  two  megaton-class  bursts  were  detonated  at  two  different 
high  altitudes,  separated  in  horizontal  distance  by  300  km,  and  in  time  by 
60  sec.  Unfortunately,  the  grid  size  outside  the  fireball  region  was  quite 
coarse  and  so  the  spatial  resolution  in  the  regions  where  molecular  emission 
might  be  important  is  very  poor.  With  this  in  mind,  we  note  that  there  is 
nothing  in  the  calculational  data  for  this  case  that  would  conflict  with 
the  conclusions  reached  in  Reference  1-1.  However,  as  emphasized  above, 
more  definitive  conclusions  must  await  the  results  anticipated  from  the 
SCENARIO  code  in  the  near  future. 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


This  report  addresses  several  topics  related  to  EXCEDE  experiments 
and,  briefly,  various  items  pertaining  to  LIVIR  emission  in  a  nuclear  environ¬ 
ment.  The  following  is  a  summary  of  the  principal  results  and  conclusions. 

Regarding  the  spatial  dimensions  and  mechanism  involved  in  the 
anomalous  region  near  the  EXCEDE : SPECTRAL  guns: 

•  A  comparison  between  the  calculated  crosswise  integrated  intensity 
and  that  inferred  from  the  ground-based  photographic  data  shows  that  the 
enhanced  optical  emission  from  the  pulse  at  123-km  altitude  occurs  in  a 
region  surrounding  the  accelerator  that  extends  for  several  hundred  meters 
along  the  magnetic  field.  A  more  accurate  definition  of  the  spatial  extent 
of  this  enhanced  region  is  not  possible,  for  the  EXCEDE: SPECTRAL  experiment, 
because  of  the  finite  sensitivity  of  the  photographic  film. 

•  An  attempt  to  salvage  the  independent  particle  model  of  energy 

o 

deposition  in  EXCEDE  experiments,  by  including  excitation  of  3914  A  radiation 
by  hot  plasma  electrons,  failed.  During  the  bombardment  time  of  a  typical 

o 

air  parcel,  the  enhancement  in  the  3914  A  emission  rate  from  this  mechanism 
is  at  most  about  50  percent.  Thus,  the  independent  particle  model  cannot 
account  for  the  observed  two  order-of-magnitude  enhancement  in  the  optical 
emission  rate.  It  appears  likely  that  collective  effects,  producing 
internal  electric  fields  that  accelerate  the  plasma  and  beam  electrons,  are 
required  to  explain  the  observations. 
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•  A  somewhat  cursory  study  of  the  beam  plasma  discharge  phenomenon 
reveals  that  the  criteria  for  its  production,  as  well  as  the  observed  optical 
consequences,  appear  to  be  fulfilled  in  the  EXCEDE: SPECTRAL  experiment. 
Although  this  does  not  provide  proof  of  the  BPD  mechanism  in  EXCEDE,  the 
circumstantial  evidence  points  strongly  in  that  direction. 

•  Application  of  the  laboratory-inferred  criteria  for  onset  of  a 
BPD  reveals  that  the  phenomenon  should  not  be  present  in  either  an  auroral 
or  a  nuclear  environment.  This  corroborates  the  assumption,  implicit  in 
auroral  and  nuclear  models  of  electron  energy  deposition,  that  collective 
effects  are  not  important.  Therefore,  the  independent  particle  assumption, 
and  the  consequences  therefrom,  should  be  essentially  valid  in  these  cases. 

Regarding  the  feasibility  of  utilizing  EXCEDE  as  a  tool  for 
verifying  our  bound-bound  models  of  IR  plasma  radiation: 

•  EXCEDE  does  not  appear  to  be  a  useful  tool  in  this  respect. 

•  If  the  excitation  and  ionization  are  produced  primarily  by 
the  beam  electrons  themselves,  then  the  bombardment  time 
required  (seconds)  for  adequate  plasma-emission  buildup 
exceeds  the  time  available,  in  a  typical  EXCEDE  experiment, 
by  one  to  two  orders  of  magnitude.  Extending  the  bombard¬ 
ment  time  by  flying  the  rocket  parallel  to  the  earth’s 
magnetic  field  is  probably  not  a  viable  option. 

•  If  BPD  conditions  exist,  the  time  requirement  should  easily 
be  fulfilled,  but  the  electron  temperature  may  remain  high 
and  severely  depress  the  bound-bound  plasma  emission  rate. 

In  any  event,  accurate  measurement  of  the  electron  tempera¬ 
ture  and  the  atomic  ion  densities  would  be  required;  measure¬ 
ments  that  are  not  likely  to  be  performed  easily  or  inexpen¬ 
sively  in  a  space  environment. 
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•  Measurement  of  the  bound-bound  plasma  radiation,  originating  from 
electron-beam  irradiation  of  dry  air  for  a  sufficient  period  of  time  (under 
non-BPD  conditions)  should  be  relatively  free  from  interference  by  competing 
molecular  emissions,  at  least  in  the  LWIR  region.  Whether  or  not  this  con¬ 
clusion  also  holds  under  BPD  conditions  is  presently  unknown. 

•  Any  electron  beam  experiment,  designed  to  investigate  bound-bound 
plasma  radiation,  is  probably  best  carried  out  in  the  laboratory. 

Regarding  some  LWIR  considerations: 

•  Excitation  of  the  infrared  bands  of  No0  by  solar  pumping  should 
be  small  compared  to  that  by  earthshine  scatter. 

•  As  anticipated,  available  data  are  insufficient  to  render  defini¬ 
tive  conclusions  regarding  the  magnitude  of  the  f  numbers  for  uranium  oxide 
ions  and,  hence,  on  the  magnitude  of  the  LWIR  emission  therefrom. 

•  Preliminary  investigations,  based  on  very  incomplete  calculational 
data,  suggest  that  the  relative  importance  of  various  LWIR  emission  processes, 
inferred  from  single-burst  calculations,  may  also  be  valid  in  a  multiburst 
environment . 
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APPENDIX  A 


ARCTIC-CODE  MODIFICATIONS  REQUIRED  FOR 
INTENSE  BOMBARDMENT  CONDITIONS 


PRELIMINARY 

A  description  of  the  ARCTIC  code,  as  written  for  use  under  auroral 
simulation  conditions,  has  been  given  in  References  2-13  through  2-16.  One 
of  the  important  quantities  computed  by  the  code  is  the  electron  temperature, 
T  .  It  is  determined  from  the  basic  equation 

dT 

sr  ’ H ' c  ■  (A-« 

where  H  and  C  are,  respectively,  the  heating  and  cooling  rates  per 
electron.  The  heating  term  arises  mainly  from  the  bombarding  flux  of 
electrons  and,  for  auroral  use,  the  cooling  terms  include  electron  energy 
loss  by:  (1)  elastic  collisions  with  N  ,  0o,  0,  and  positive  ions,  (2) 

3 

rotational  excitation  of  N?  and  CL,  (3)  fine  structure  excitation  of  0('P), 

Z  Z  11 
(4)  vibrational  excitation  of  N2  and  0o,  (5)  excitation  of  0(  D)  and  0(  S) 

from  0(^P),  (6)  excitation  of  0+(^D)  and  0+(^P)  from  0+(^S),  and  (7)  excita- 
2  2  4 

tion  of  N(  D)  and  N(  P)  from  N(  S)  .  Deexcitation  of  the  above  excited  states 
by  plasma  electrons  gives  a  negative  contribution  to  the  cooling  rate  which, 
in  effect,  adds  to  the  electron  heating  rate. 

The  foregoing  cooling  mechanisms  are  entirely  adequate  for  auroral 
conditions  in  which  the  electron  temperature  seldom  exceeds  about  3000  °K. 
However,  under  the  intense  bombardment  conditions  encountered  in  an  EXCEDE 
experiment,  the  electron  temperature  can  become  sufficiently  high,  at  least 
at  early  times  after  beam  turn  on,  when  the  plasma  density  is  still  quite  low, 
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to  excite  some  of  the  higher-lying  electronic  states  of  the  air  atoms, 
molecules,  and  ions.  Consequently,  for  an  adequate  representation  of  the 
electron  temperature  under  such  conditions,  as  well  as  for  a  determination 

o 

of  the  (plasma)  electron-impact  excitation  rate  of  3914  A  emission,  it  is 
necessary  to  include  these  other  states  in  the  cooling  process. 


Furthermore,  a  review  of  the  models  in  ARCTIC,  valid  for  auroral 
applications,  revealed  that  those  for  electron  cooling  by  elastic  collisions 
with  Nn,  0o,  0,  and  by  vibrational  excitation  of  N0  and  0o,  are  not  valid 
for  large  values  of  T  .  Consequently,  it  was  necessary  to  revise  them. 

This  appendix  describes  the  additional  cooling  processes  required, 

and  presents  analytical  expressions  for  their  rates  together  with  revised 

rates  for  the  other  processes  mentioned  above.  Also  presented  are  the  rates 
0 

for  3914  A  emission  and  for  ionization  produced  by  hot  plasma  electrons. 


REVISION  OF  ELECTRON  COOLING  RATES 


Elastic  Scattering  by  N^,  0^,  0 

The  cooling  rate  per  electron  due  to  elastic  collisions  with  gas 
molecules  can  be  written,  in  units  of  ergs  electron  ^  sec  \  as  (Reference 


A- 1)  : 


I  dE\ 

{dt) 


=  -  4N 


elastic 


mM  ,  / 8k  \1//-/Te  ,  (m)^ 


k(T  -T) 


(m+M) 


2  v  e 


n ■ 


m 


+  M  / 


(A-2) 


where 


/a(m) 


o 

>■>■/ 


5  (m)  ,  .  -Kv  . 
vo  (vj  e  dv 


(A-3) 


and 


(A-4) 
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Here,  m  and  M  are  the  masses  of  the  electron  and  gas  molecule,  respec¬ 
tively,  T  is  the  neutral  gas  temperature,  and  a^(v)  is  the  momentum 
transfer  cross  section  for  electrons  of  velocity  v  in  a  gas  of  number 
density  N. 


expression 


where  E 


We  have  approximated  the  momentum  transfer  cross  sections  with  an 
of  the  form 


m ^1/2.  -bE  ,  2. 

a  (E)  =  0Q(l  +  aE  )e  (cm  ) 


(A-5) 


is  the  electron  energy  (eV) .  Alternatively, 


(m) ,  ,  in-8  ,  -2.8x10  16v2 

a  (v)  =  0^(1+1.69*10  av)e 


(A-6) 


Equations  A- 2  through  A-6,  together  with  the  approximation 


mM 

(m+M) 2 


m  _  5. 48  x  10 

M  .  4< 


where t U  is  the  molecular  weight  of  the  gas,  lead  to  the  cooling  rate 
(eV  electron  ^  sec  b : 


elastic 


0  N  1  /?  r  - 

0. 118  T  '  (T  -T) (1+8.618X10  bT  ) 
e  e  e 


1.54  x 10 


l+8.6lxio"5hT 


(A-7) 


The  cross  sections.  Equation  A-5,  have  been  fitted  to  data  for  N0 , 
0^,  and  0;  the  fits  being  good  for  plasma  electron  energies  up  to  about 
10  eV.  Values  for  the  cross  section  parameters  a  ,  a,  and  b,  along  with 
references  to  the  data,  are  shown  in  Table  A-l. 
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Table  A-l.  Parameters  used  in  Equation  A-5  to  fit  measured  elastic 
scattering  cross  sections  for  Ng,  0^  and  0. 


Mol ecul e 

a 

0 

(an^) 

a 

lev-’/2) 

b 

(eV 1) 

References 

N2 

28 

X 

O 

1 

cn 

10.39 

3. 38  x  10"2 

A-2 ,  A- 3,  A-4 

°2 

32 

3x  10’ 17 

3.54 

1 . 68  x  10"2 

A- 3,  A-5 

0 

16 

7x  10"17 

4.0 

2.57x  10"2 

A- 3,  A- 6 

Vibrational  Excitation/Deexcitation  of 

Involved  here  is  the  collection  of  reactions  of  the  form 


e  +  (v) — ►e  +  N^fv')  ;  v'  >  v  (A-8) 

with  forward  and  reverse  rate  coefficients  k  ,  and  k  ,  ,  respectively. 

vv  V  V 

Thus,  the  net  cooling  rate  per  electron  is 


(-) 

WL 


N„(vib.)  v,v 


£,  |kv'vCN,2(v,;i]  ‘  [No(v)]|w,rt 


vv*  2  ^  J  J  |  v *  v 


(A-9) 


where  W  ^  is  the  energy  difference  between  vibrational  states  v*  and  v. 


Equation  A- 9  has  always  been  an  integral  part  of  the  electron 
temperature  model  in  the  ARCTIC  code.  However,  the  simplified  rate  coef¬ 
ficients  used  for  auroral  applications  were  found  to  be  inapplicable  at 
the  high  values  for  T^  encountered  in  EXCEDE  experiments.  Extensive 
modifications  have,  therefore,  been  made  using  the  measured  cross  sections 
of  Schultz  (Reference  A- 7)  and  calculations  by  Chen  (Reference  A-8). 
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Since  the  vibrational  excitation  of  for  electron  energies 
^1.6  eV  proceeds  through  formation  of  temporary  negative  ion  states,  each 
cross  section  consists  of  a  series  of  resonance  peaks.  In  order  to  facili 
tate  the  analytical  representation  of  each  reaction  rate  coefficient,  we 
have  approximated  the  resonance  peaks  by  a  series  of  sawteeth.  The  rate 
coefficient  for  each  forward  reaction  can  then  be  expressed  in  the  form 


7  2x10'  {  2  a(E  +AE  )  a(E  -AE  ) 

/.2  x  1U  ^  i_  )  rn  xAp  3  J  X  tv  .Ar  'i2^  J  3 


V,  V  ' 


T  T  ji-  J  (E.+AE.fe  3  3  +  (E.-AE.)‘ 

T  1/2  *—>  ae.  )  *■  ]  r  1  r 

e  1  J  ( 

aE. 


2E 


5.  e  3  +  -=•  (2+aE. 


J 


(2  +  aEj  )  (aF.j  - 1)  e 


aE. 
1 


a (e . +AE . ) 

[2+a(E.+AE.)l [a(E.+AE.)-l]e  3  3 

I  J  J  l  


a  (E .  -AE .  ) 

[2+a(E.-AE.)l [a(E.-AE.)-l]e  3  3 

J  3  1  J 


(A- 10) 


where 


a  =  -  1 . 16  x  10  /T 


(A- 11) 


and  the  summation  is  over  the  resonances  for  each  of  the  excitation  cross 


sections.  Values  of  the  parameters  a.,  E.,  and  AE.  for  each  excitation 
reaction  are  shown  in  Table  A-2. 


Values  for  the  inverse  (deexcitation)  rate  coefficients  can  be 
obtained  (using  detailed  balancing)  from  the  equation 


v'  v 


=  e 


W  ,  /kT 
vf  v  e 


vv' 


V!  >  V 


(A-12) 


From  the  known  vibrational  energy  levels  of  N0  (Reference  A-9) ,  we  can  sub 

stitute  in  Equation  A-12  for  W  .  and  obtain  the  result 
n  v?  v 

kv,v  =  exp{3.38xl03(v'-v)[l-6.13xl0"3(l+v+v')]/Te}kvv,  (A-13) 
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Table  A-2.  Values  of  parameters  for  use  in  Equation  A- 1 0 


k*v’ 

rl 

El 

^E1 

'2 

“2 

iE2 

*3 

c 

-3 

1E3 

J4 

E4 

iE4 

°5 

E5 

iE5 

^6 

E6 

iE6 

k0  -  1 

0.38 

1.5 

1.1 

1,5 

2.03 

0.2 

1.4 

2.32 

0.18 

1.05 

2.62 

0.19 

0.5 

2.83 

0.12 

0.4 

3.12 

0.29 

’*0  -  2 

1.1 

2.08 

0.25 

1.0 

2.44 

0.18 

0.56 

2.75 

0.15 

0.32 

3.0 

0.18 

0.12 

3.31 

0.38 

0 

0 

1 

k0-3 

0.9 

2.17 

0.29 

0.5 

2.61 

0.20 

0.30 

2.92 

0.17 

0.12 

3.21 

0.15 

0 

0 

1 

0 

0 

1 

*0  -  4 

0.8 

2.25 

0.29 

0.4 

2.72 

0.15 

0.16 

3.04 

0.19 

0 

0 

1 

0 

0 

1 

0 

0 

1 

k0  -  5 

0.56 

2.37 

0.31 

0.21 

2. 35 

0.19 

0.14 

3.19 

0.18 

0 

0 

1 

0 

0 

1 

0 

0 

1 

k0  -  6 

0.51 

2.54 

0.3 

0.17 

3.00 

0.22 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

k0  -  ? 

0.32 

2.69 

0.3 

0.09 

3.19 

0.23 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

kl  *  2 

0.25 

1.54 

0.6 

3.0 

1.71 

0.2 

0,3 

2.46 

0.2 

0.3 

2.78 

0.25 

0.25 

3.0 

0.18 

0.18 

3.23 

0.2 

kl  -  3 

1.3 

1.8 

0.3 

0.16 

2.37 

0.25 

0.3S 

2.61 

0. 16 

0.24 

2.84 

0.16 

0.10 

3.2 

0.25 

0.08 

3.4 

0.1 

kl  *  4 

0.5 

1.86 

0.3 

0.2 

2.4 

0.25 

0.21 

2.6 

0.22 

0.16 

2. 35 

0.15 

0.09 

3.2 

0.14 

0 

0 

1 

k1  *  5 

0.2 

’  .92 

0.3 

0.07 

2.38 

0  4 

0.25 

2.57 

0.2 

0.18 

2.35 

0.16 

0.08 

3.2 

0.22 

0 

0 

1 

k!  -  6 

0.05 

1.35 

0.15 

0.15 

2.36 

0.23 

0.16 

2.63 

0.12 

0.15 

2.78 

0.12 

0.15 

3.0 

0.15 

0.08 

3.3 

0.15 

kl  -  7 

0.2 

2.41 

0.18 

0.12 

2.59 

0.11 

0.16 

2.33 

0.18 

0.08 

3.2 

0.20 

0.05 

3.45 

0.08 

0 

0 

1 

k2  -  3 

0.15 

1.26 

0.35 

1  6 

1  42 

0.17 

0.4 

1.58 

0.19 

0.5 

2.0 

0.22 

0.2 

3.03 

0.8 

0 

0 

1 

*2-4 

0.4 

1.41 

0.2 

0.6 

1,78 

0.28 

0.25 

2.1 

0.18 

0.22 

2.9 

0.6 

0 

0 

1 

0 

0 

1 

k2  -  5 

0.8 

1.35 

0.2 

0.1 

2.18 

0.3 

0.18 

3.0 

0.6 

0 

0 

i 

0 

0 

1 

0 

0 

1 

k2  -  6 

0.54 

1.91 

0.23 

0.15 

2.76 

0.12 

0.  14 

3.0 

0.11 

0.09 

3.3 

0. 14 

0.08 

3.5 

0.1 

0 

0 

1 

k2  -  7 

0.23 

1.92 

0.15 

0,23 

2.08 

0.15 

0.15 

3.05 

0.6 

0 

0 

1 

0 

0 

1 

0 

0 

1 

k3  -  4 

0.12 

1.0 

0.5 

0.55 

1.12 

0.23 

2.0 

1.4 

0.20 

0.3 

1.87 

0.19 

0.2 

2.2 

0.25 

0 

0 

1 

k3  -  5 

1.0 

1  46 

0.25 

0.4 

2.3 

0  35 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

*3  -  6 

0.29 

1.42 

0.1 

0.29 

1.6 

0.15 

0.15 

1.82 

0. '  6 

0.25 

2.08 

0.23 

0.15 

3.26 

0.41 

0 

0 

1 

k3  -  7 

0.1 

1.58 

0.23 

0.22 

1.85 

0.15 

0.12 

2.04 

0.11 

0.12 

2.25 

0.15 

0.08 

2  89 

0.08 

0.09 

3,38 

0.53 

k4  -  5 

3.2 

0  95 

0.5 

1.5 

1.06 

0.15 

0.7 

1.31 

0.13 

0.18 

1.57 

0.35 

0.18 

2.02 

0.15 

0.11 

2.35 

0.17 

k4  -  5 

0.4 

1.03 

0.17 

0.71 

1  28 

0.15 

0  69 

1.43 

O.H 

0.15 

1 . 77 

0.15 

0.31 

2.1 

0.  ’6 

0 

0 

i 

k4  -  7 

0.45 

1.25 

0.17 

0.74 

1.46 

0.17 

0.19 

1.92 

0.15 

0.15 

2.22 

0.16 

0 

0 

1 

0 

0 

1 
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Vibrational  Excitation  of  0^ 


Using  cross  sections  measured  by  Hake  and  Phelps  (Reference  A-5), 
we  have  computed  the  electron  cooling  rate  (eV  electron  *  sec  b>  valid  to 
high  values  of  T  .  The  result  is  as  follows: 


dE\ 

'C>2  (vib. ) 


10"5[0„1  /  -4700/T  -6875/T  -9050/T 

2  /  e  0  e  0  _  e 

T  3/2  >e  *  2e  *  2-30 

e 

-11210/T  -13400/T  - 15570/T 

7. 9e  e  +  15. 5e  e  +  15. 9e 


-17700/T  -20300/T 

+  14. le  e  +  26. le 


(A- 14) 


The  eight  terms  correspond,  respectively,  to  contributions  from  excitation 
of  the  first  eight  vibrational  levels 

COOLING  BY  ELECTRONIC  EXCITATION 
AND  IONIZATION  OF  N2,  02,  0 

Table  A- 3  shows  the  additional  electronic  states  of  09,  0, 

and  their  ions,  along  with  their  energies  relative  to  the  ground  state  of 

the  neutral  molecule,  that  we  now  consider  as  candidates  for  excitation  by 
* 

plasma  electrons.  This  subsection  presents  the  analytical  formulas  that 
have  now  been  incorporated  into  ARCTIC  to  calculate  the  plasma-electron 
cooling  rates  by  excitation  of  these  additional  states. 


Let  k. .  be  the  rate  coefficient  for  electron  excitation  of 
ij 

state  j,  at  energy  W.,  from  the  ground  state  of  species  i.  Then  the 


*  These  states,  among  many  others,  have  always  been  an  integral  part  of 

the  ARCTIC  code  insofar  as  their  excitation  by  beam  electrons  is  concerned. 
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Table  A-3.  Additional  electronic  states  included  in 
electron  temperature  calculation. 


State 

Energy  (eV) 

State 

Energy  { eV ) 

State 

Energy  (eV) 

(1)  N2(aV) 

6.2 

(9)  N2+(D2ng) 

22.0 

(17)  o2+(A2nu) 

16.9 

(2)  N2(B3ng) 

7.4 

(10)  N2+(cV) 

23.6 

(18)  02+(b4r") 

18.2 

(3)  N2(c3nu) 

n.o 

(11)  02(a’Ag) 

1.0 

(19)  0^+ ( compos i te) 

23.0 

(4)  N2(a1Ilg) 

8.6 

(12)  o2(b1i:g) 

1.64 

(20)  0{ Rydberg) 

9.2 

(5)  N2(b1nu) 

12.8 

(13)  02(aV) 

4.5 

(21)  O(Rydherg) 

9.6 

(6)  N2+(xV) 

15.6 

(14)  02(83E”) 

8.4 

(22)  0+(4S) 

13.6 

(7)  N2+(A2nu) 

16.7 

(15)  o2+(x2ng) 

12.1 

(23)  0+(?D) 

16.9 

(8)  N^(bV) 

18.75 

(16)  02+(a\) 

16.1 

(24)  0f(2P) 

18.5 

cooling  rate,  per  electron,  due  to  excitation  of  species  i,  is 


l~)  =  -  [i]  £k..W. 

\dt  /.  1  J  ^  in 

'  !l  3 


(A- 15) 


where  [i]  is  the  concentration  of  species  i. 
given  by 


But  the  rate  coefficient  is 


kij  =  (Vaj(i:i)  =  Vi  fEU2(J j(ihE)f(E)dE  ,  (A- 16) 

W. 

J 

where  a  ^(E)  is  the  cross  section  for  excitation  of  state  j  of  species 
i  by  an  electron  of  energy  E,  and  f(E)  is  the  electron  energy  distribution 
function . 


For  a  Boltzmann  distribution  of  electron  energies. 


f(E)  = 


2tt 


(wkTJ 


3/2 


E1/2e-E/kT 


(A- 17) 
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and  Equation  A- 16  can  then  be  written,  with  E  expressed  in 

oo 


k.  . 


dE 


eV  units,  as 
(A-18) 


4 

where  the  parameter  b  has  the  value  1 . 159 x  10  .  Equations  A-15  and  A-18 
can  now  be  combined  to  yield  the  following  equation  for  the  cooling  rate 
(eV  electron  ^  sec  by  electronic  excitation  of  species  i: 


8. 36  x  1 013  [i ] 

_  3/2 


j  w. 


(E)  Ee 


-bE/T 


dE 


(A-  1 9) 


The  cross  sections  employed  in  the  ARCTIC  code  (Reference  2-13) 
for  excitation  of  the  states  shown  in  Table  A- 3  have  been  used  in  Equation 
A-19  to  obtain  cooling  rates  for  electronic  excitation  of  N?>  0?,  and  0.  In¬ 
cluded  are  states  1  through  10  in  Table  A-3  for  N?,  states  11  through  19  for  0o, 
and  states  20  through  24  for  0.  The  integrals  were  then  performed  numerically, 
for  selected  values  of  T  ,  and  the  results  fitted  with  the  following  analy- 

. e  -i  - l 

tical  expressions,  each  in  units  of  eV  electron  sec 


=  -3.92  x  10"°[N.Jexp{-  [6.  04xl04/Te+0.  16  (JlnT  -  15 . 42)  2]  } 

(A-20) 

=  -2. 45  x 10'6[02]exp{-[1.662xl03/T  0'6578 

0o (el ectronic) 

+  0. 1528(£nTe-15.42)2]}  (A-21) 

2.  86  x  10"6[0]exp{-[2.437xl05/Te1,086+0.  17  (SjiT  - 15 . 42)2]  }  . 

(A-22) 


(-) 

04  ^O(electronic) 


f  dE  \ 
idt ) 


N_ (electronic) 
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IONIZATION  BY  HOT  PLASMA  ELECTRONS 


Excitation  of  states  6  through  10  in  Table  A-3  by  plasma-electron 
impact  on  N0(X)  produces  additional  ionization  of  N0.  Similarly,  excitation 
of  states  15  through  19  and  states  22  through  24  produces  additional  ioniza¬ 
tion  of  0^  and  0,  respectively.  For  example,  the  ionization  rate  of  ^  by 
the  process 


e  +  (X)  -*■  2e  + 


(A- 23) 


can  be  written  as 
dN, 


-ar-  kN;W 


(A-24) 


where  k^T  +  is  the  rate  coefficient  for  the  process. 

Equation  A- 18  has  been  used,  together  with  the  cross  sections  in 

the  ARCTIC  code  for  excitation  of  states  6  through  10,  to  obtain  values  for 

k.T  ,  at  selected  values  of  T  .  The  result  has  been  fitted  analytically 
N2  e  3  -i 

and  can  be  written,  in  units  of  cm'  sec  ,  as 

k  +  =  2.06x  lo"7exp{- [8.436xl05/T  1  * 167+0. 16(£nT  -15. 42)2]}  (A-25) 

N2  e  e 

Corresponding  expressions  for  the  ionization  rate  coefficients  of  and  0, 
obtained  in  a  similar  manner,  are  as  follows: 

k  +  =  1. 35  x  l0_7exp{- [2. 643xl05/T  1  ‘  059+0. 1528 (S,nT  -15. 42)2]}  (A-26) 
°2  e  e 

k  +  =  1.76x  io'?exp{-[1.2xl06/T  1,225+0.17(£nT  -15.42)2]}  (A-27) 

0  e  e 


Equations  A-25  through  A-27  have  been  incorporated  into  the  ARCTIC 
code  and  provide  additional  ionization  whenever  Te  is  sufficiently  high. 
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PLASMA-EXCITATION  OF  N2+  FIRST  NEGATIVE  BANDS 


In  addition  to  excitation  of  the  First  Negative  band  system 
by  a  primary  electron  flux,  excitation  by  plasma  electrons  may  also  con¬ 
tribute  provided  is  high  enough.  In  particular,  excitation  may  occur 

by  impact  on  ground-state  molecules: 


e  +  N-fxV)  +  2e  +  N„+(B2E+) 

2K  gJ  2  K  uJ 

or  by  impact  on  ground-state  N0+  ions: 

e  +  N„+(X2E+)  -*■  e  +  N„+CB2Z+) 

2  g  2  u 


(A-28) 


(A-29) 


Equation  A-18  has  been  used,  with  the  ARCTIC-code  cross  section 

for  process  A-28,  to  compute  values  for  the  rate  coefficient,  k00,  at 

28 

selected  values  of  T  .  The  result  has  been  fitted  analytically  and  can 

e  3  -1 

be  written,  in  units  of  cm  sec  ,  as 


k  =  2.55 x  l0"8exp{- [4.83xl05/T  1,095+0. 157(£nT\-15.42)2]}  .  (A-30) 

28  6  e 


Similarly,  using  the  cross  section  for  process  A-29  measured  by  Crandall, 
et  al.  (Reference  A-10),  we  find 


k  =  2.66  x  io'8exp{-[7.1xl06/Te1’757+0.3408(£nTe-11.513)2]}  .  (A-31) 


It  can  be  shown  (Reference  2-13)  that  the  volume  emission  rate 

°  + 
of  3914  A  photons  from  the  (0-0)  band  of  the  N  First  Negative  system  is 

given  by 

l>3914  «  0.65  P[N2+(B)]  (A-32) 
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where  P [N2  (B)]  is  the  production  rate  of  the  state  of  N^+.  Conse¬ 

quently,  the  volume  emission  rate  of  3914  A  photons,  excited  by  hot  plasma 
electrons,  is 


(plasma) 

$ 

3914 


°.65  Ne(k28[N2]+k29[N2 


]) 


(A- 33) 


Equations  A-30  through  A-33  have  been  incorporated  into  the 
ARCTIC  code. 
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